arXiv:1507.07935vl [astro-ph.GA] 28Jul2015 


Draft version July 30, 2015 

Preprint typeset using D-TfTjX style emulateapj v. 5/2/11 


RAM PRESSURE STRIPPING OF THE LARGE MAGELLANIC CLOUD’S DISK AS A PROBE OF THE 

MILKY WAY’S CIRCUMGALACTIC MEDIUM 

Munier Salem 1 , Gurtina Besla 2 , Greg Bryan 1 , Mary Putman 1 , Roeland P. van der Marel 3 , Stephanie 

Tonnesen 4 

Draft version July 30, 2015 

ABSTRACT 

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), 
implying a well-constrained pericentric passage about the Milky Way (MW) ~ 50 Myr ago. In 
this scenario, the LMC’s gaseous disk has recently experienced stripping, suggesting the current 
extent of its HI disk directly probes the medium in which it is moving. From the observed stellar 
and HI distributions of the system we find evidence of a truncated gas profile along the windward 
“leading edge’ of the LMC disk, despite a far more extended stellar component. We explore the 
implications of this ram pressure stripping signature, using both analytic prescriptions and full 3- 
dimensional hydrodynamic simulations of the LMC. Our simulations subject the system to a headwind 
whose velocity components correspond directly to the recent orbital history of the LMC. We vary the 
density of this headwind, using a variety of sampled parameters for a /3-profile for a theoretical 
MW circumgalactic medium (CGM), comparing the resulting HI morphology directly to observations 
of the LMC HI and stellar components. This model can match the radial extent of the LMC’s 
leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is 
n p (R = 48.2 ± 5 kpc) = l.V'45 x 10 -4 cm -3 . The implied pericentric density proves insensitive 
to both the broader CGM structure and temperature profile, thus providing a mo del-independent 
constraint on the local gas density. This result imposes an important constraint on the density profile 
of the MW’s CGM, and thus the total baryon content of the MW. From our work, assuming a /3-profile 
valid to ~ r v j r , we infer a total diffuse CGM mass M(300 kpc) = 2.6 ± 1.4 x 10 10 Mq or approximately 
15% of a 10 12 M q MW’s baryonic mass budget. 

Subject headings: (galaxies:) Magellanic Clouds — Galaxy: structure — galaxies: kinematics and 
dynamics - hydrodynamics — (galaxies:) intergalactic medium 


1. INTRODUCTION 

The Circumgalactic Medium (CGM) of the Milky 
Way (MW) represents a gaseous, multiphase plasma 
surrounding the Galactic plane and extending out to 
the virial radius. Though substantially more rar¬ 
efied than the interstellar-medium (ISM), the density 
of plasma within this vast volume remains relatively 
uncertain. Studies in X-ray have detected CGM ma¬ 
terial ~ 10 6 K, though the total mass, and whether 
or not it accounts for a majority of baryons in the 


MW 

remains in dispute dWans: & McCrav 

199.1 1 Wan g 

2005 

IGuDta et all 2012t Mathur 20121 

Wang & Yao 

2012 
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of L ~ A* galaxies at low-redshift have likewise found 
evidence for a large fraction of galactic baryons residin g 
in the CGM (jTumlinson et alJl201lHPeeples et al.ll2014h . 
a majority of which may be in the for m of diffuse gas at 
surprisingly low tempe ratures < 10 5 K dWerk et al.ll20l4 
iTumlinson et al.ll2013H . 

The CGM represents an important store of fuel for 
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star formation (SF) as well as a tracer of inward flows 
from the pristine intergalactic medium (IGM) and out¬ 
ward flows of e nriched material from the star-forming 
disk plane ( e.g. lBinnevlll977t Ivan de Voort et al.l 1201 lfc 


Cowie et al. 1 19951: Dave et all 120121 iFraternali et all 
20131 ISuresh et aLl 12015T) . Thus characterizing its com¬ 


position aids in a comprehensive understanding of galaxy 
formation. For a recent review of the many components 
of the CGM, see lPutman et al.l (12012H . 

Hydrodynamic probes of the MW CGM gas structure 
provide a powerful compliment to observational studies 
in absorption and emission. Modeling ionization states 
relies on assumptions regarding the gas’ temperature and 
sources of photoionization; a single waveband study can¬ 
not probe the broad range of gas phases spanned by the 
CGM in all its forms. In contrast, studies that exploit 
ram pressure stripping (RPS) probe diffuse gas of all 
temperatures and ionization states, and the stripping dy¬ 
namics a re insensitive to the oncom ing wind’s tempera¬ 
ture fe.g. lRoediger fc Hensleril2005L find only a slight de- 
pe ndence on the Mach n u mber of the flow). Autho r s such 
as [Grcevich fc Putniaiil (120091) and iGatto et al.l (j2013l ) 
considered a host of dwarf spheroidals orbiting the MW, 
most devoid of HI, to place preliminary bounds on the 
CGM’s gas profile. In this paper we consider instead 
the Large Magellanic Cloud (LMC), a relatively mas¬ 
sive, late-type dwarf, with well-studied HI and stellar 
disks and well-constrained orbital properties. 

The Magellanic Clouds (MCs) move through the MW 
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at ~ 50 kpc from the Galactic Center at velocities ~ 300 
km/s, with a Leading Arm (LA) and trailing Magel¬ 
lanic Stream (MS) of gas strewn across much of the 
southern sky. The interaction of this complex with the 
CGM likely alters the appearance and dynamics of these 
gaseous components. The LMC’s proper motion vector 
implies that any CGM headwind would most directly 
impact the disk’s northeastern edge, referred to here¬ 
after as its “leading edge”. The HI profile here truncates 
abruptly at R ~ 6.2 kpc from the kinematic center of the 
LMC (see Sec tion 12. II) . in contrast to other quadrants of 
the LMC ('e.g. lPutman et afll2003D . Despite the absence 
of gas, the stellar profile continues uninterrupt ed well 
beyond this radius (jvan der Marel fc Cionill2001l) . which 
would rule out a tidal explanation for this HI truncation 
(as tides would truncate both gas and stars). The lead¬ 
ing edge is characterized by multiple HI velocity com¬ 
ponents , the faster among t hem being possibly extra- 
planar dLuks fc Rohlfsl IT992t iNidever et al .1120081I 2010 D . 
which point towards a strong interaction between the 
LMC’s ISM and the ambient material it passes through. 
HI column and star formation are both more concen¬ 
trated in the southeastern portion of the LMC disk; an 
entire “supershell” of denser gas and rapid star forma¬ 
tion exists at the lea ding edge of the LMC disk, perhap s 
due to a bow-shock (Ide Boer et al.|[l998t lMuralill2000i ), 
and regions devoid of gas exist despite t he presence of 
young stars (llndu fc Subramaniamll 20 Tll ). Likewise the 
MS and LA are likely significantly altered by this am- 
bient medium. The stream exhibits strong Ha emission 
dBland-Hawthorn et al.l 120071 : IWeiner fc Williams! Il996l) 
and high velocity clouds (HVCs) kinematically associ¬ 
ated with the whole complex consistently feature head- 
tail structures well aligned with t he orbital movement 
of the Magellan ic Complex at large (|Putman et al.ll 20 TTl : 
iFor et aUPOTl . 

Precise proper motion (PM) measurements of MC 
stars in the past decade have substantially im proved our 
know l edge of the LMC’s PM vector dKallivavalil et all 
l200f& iPiatek et al.1 120081 1 and enhanced our under¬ 
standing of the LMC’s orientation and internal rota- 
tion and the dynamical re lationship between the MCs 
dKallivavalil et akl 120061 : Ivan der Marel fc Kallivavalill 
120141) . These results imply that the LMC’s kinetic en¬ 
ergy is on par with its binding energy to the MW, leading 
to a picture in which the LMC passed through its last 
perigalacticon less than 100 Myr ago and is either on its 
first infall or else on a highly eccentric orbit whose pre¬ 
vious perigalactic on occurred a sub stantial fraction of a 
Hubble time ago (|Besla et al.ll2007l hereafter B07). The 
most recent epoch of PM measurements slightly reduced 
the implied LMC proper motion velocity, but still con- 
firms this h ighly eccent ric, possibly first infall scenario 
(| K alii vavalil et, all 1 2 01 3l hereafter K13). 

A detailed study of ram pressure stripping’s effect 
on the LMC disk gas’ e xtent was last undertaken by 
iMastropietro et ahl ( 20051) who followed the LMC’s or¬ 
bital history within an N-body/SPH simulation of both 
the MW and LMC’s dark and baryonic components. 
Their work suggested a MW gaseous halo density of 
« 8 x 10 -4 cm ' 3 at 50 kpc from the Galactic center. 
However, they used a substantially lower-energy LMC 
orbit, which allowed for multiple pericentric passages of 
comparable depth into the MW potential, and a much 


longer interaction time between the LMC gaseous disk 
and the MW CGM, during which time the headwind 
compressed and heated the LMC disk material, allow¬ 
ing for slow, continual stripping down to the present ob¬ 
served truncation radius — a decidedly delocalized probe 
of this ambient material. The slower orbit also enhanced 
the role of tidal interactio ns between the MW and LMC. 
IMastropietro et akl (120091) used the updated orbital sce¬ 
nario of B07 but focused on the effects of compression- 
induced star formation and did not vary the halo density 
or consider h ow the leading edg e infers a MW diffuse 
CGM density. iBesla et ah (2012|) explored a toy simula¬ 
tion of the evolution of the gaseous Magellanic System 
wherein a ram pressure acceleration term is applied uni- 
formaly to all gas particles. They find that the struc¬ 
ture of the LMC’s gas disk should be very sensitive to 
the properties of the ambient medium through which the 
galaxy is moving. 

The new, longer period, eccentric orbits of B07 and 
K13 no longer allow for slow, continual stripping mech¬ 
anisms, and suggest a scenario in which the current 
HI truncation radius along the leading edge may be 
set by fast st r ipping , well described by the model of 
I Gunn fc Gott (jl972f ) (hereafter G72), having occurred 
during the LMC’s recent pericentric passage of unprece¬ 
dented speed and depth within the MW potential. This 
points to a picture where the LMC’s current HI extent 
along the leading edge may provide a direct probe of the 
MW CGM density, likely localized to simply a measure¬ 
ment of this density at the LMC’s recent perigalacticon 
at R = 48.2 ± 2.5 kpc from the Galactic Center. Corrob¬ 
orating this notion and deriving a bound for this density 
is the primary focus of this paper. 

This paper is organized as follows. In Section [2] we 
introduce an analytic model for the LMC internal struc¬ 
ture and dynamics, the LMC’s motion through the MW 
halo and the structure of the MW CGM. In Section [3] 
we use this setup along with an analytic description of 
RPS to construct a toy model through which we arrive 
at a localized density constraint at r ~ 50 kpc. In Sec¬ 
tion [4] we then explore this process in greater detail via 
three-dimensional hydrodynamic simulations of the LMC 
passing through the MW halo. Section [5] synthesizes the 
analytic and simulated results together along with an 
enumeration of relevant uncertainties to form an accu¬ 
rate diffuse gas density bound. Section [ 6 ] discusses the 
implications of our work for the total CGM mass bud¬ 
get and how our work informs the current understanding 
of the Magellanic Stream’s (MS’s) formation. Section [7] 
summarizes our main findings. 

2. THEORETICAL MODEL FOR RPS OF THE LMC 

In this section, we outline a model framework that 
will be applied consistently to all analytic and simulation 
work in this paper. Sections 12.11 12.31 outline our pre¬ 

scriptions for the LMC internal structure and mass com¬ 
ponents, the MW diffuse gaseous halo, and the LMC’s 
center of mass (COM) motion through the MW, respec¬ 
tively. 


2 . 1 . Model for the LMC 

We begin by describing our model for the LMC’s bary¬ 
onic and DM distributions. Table 1 summarizes the im¬ 
portant parameters. 
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Static DM Potential 

Spherical Profile (Burkert 1995) 
po 3.4 x 10 -24 g/cm 3 

ro 3 kpc 

M(100 kpc) 5 x 1O 1O M 0 

Static Stellar Potential 

Plummer-Kuzmin Disk 
- Mi^arnoto - ^_^Jagai (1975) 

M* 2.7 x 10 9 Mq 

a* 1.7 kpc 

£>* .34 kpc 

Initial Gas Distribution 

Exponential Disk 
- Tonnesen_&_^r^an (2009) 

M gas 5 x 10 s Mq 

a gas 1.7 kpc 

6 gas -34 kpc 

Mrot 7.2 x 10 8 Mq 


Dynamical Mass 

M dyn (8.7 kpc) 1.4 x 1O 1O M 0 


Table 1: Parameters for our LMC model, which consists of 
static DM and stellar potentials and a self-gravitating gas disk 
with an initially exponential profile. For this simple setup, 
all profiles share a common kinematic center and angular mo¬ 
mentum axis, and the orbits are completely circularized. Here 
the dynamical mass includes all components (stars, gas and 
DM). 


For the stellar an d gaseous disks of t h e LM C we 
follow the setup o f IRoediger fc Briiggeul (120061) and 
iTonnesen fc Brvanl (i2009h . hereafter T09. We assume 
the stellar and gaseous disks are circularized, share a 
common kinematic center, and are aligned in their rota¬ 
tion axes. In reality, the present day stellar distribution 
is intrinsically asymmetric and elliptical, and the position 
angle of progressively fain ter isop h otes of the stel l ar dis - 
tribution rotate by ~ 90° (|van der Marel fc Cionil[200lD . 
However, these simplifications should not substantially 
alter our RPS dynamics, which is mostly governed by the 
stellar and gaseous surface densities at R ~ 6 kpc and to 
a lesser extent the circular velocity of the HI , related t o 
the total enclosed mass at these radii. iKim et abl ( 1998 ). 
hereafter K98, studied the HI velocity field to find an 
upper limit on the dynamical mass within a radius of 4 
kpc of M (4 kpc) = 3.5 x 1O 9 M 0 . From their HI surface 
maps, they also inferred a total disk mass of 2.5 x 1O 9 M 0 
and a gaseous component of 5 x 1O 8 M 0 within the same 
radius. We use these values to constrain the normaliza¬ 
tion of our profiles. F or the stars we emp l oy a static 
Plummer-Kuzmin disk (iMivamoto fc Nagailll975l) . 


$*(R, z) = GM+ 



(a* + y/z^+b f) 



(1) 


with mass M* = 2.7 x 1O 9 M 0 , radial scale length a* = 
1.7 kpc and vertical scale height b+ = .34 kpc. T hese 
scale radii are motivated by Ivan der Marel fc Cionil 
(j2001l l. For the gas disk, we adopt the profile of T09, 


3 {R,z) = ^ sah —0.5 2 sech 

27 ™gaAas 


R \ 


j sech 

^gas / 


gas 


( 2 ) 


where a gas and & gas are the radial and vertical scale 


heights of the disk and M gas is a normalization related 
to the total mass of the disk by Mp, gas ~ 1.44M gas . We 
consider only profiles that roughly obey the constraint 
of K98: M gas (4 kpc) = 5 x 1O 8 M 0 . For the majority of 
simulations, we then set the gaseous scale radii equal to 
that of the stellar distribution, and thus a gas = 1.7 kpc, 
bgas = -34 kpc and M gas = 5 x 1O 8 M 0 . This “charac¬ 
teristic mass” is identical to the mass within 4 kpc by 
coincidence; the total enclosed mass of this disk model 
to large radii Mt ~ 1.44M gas = 7.2 x 1O 8 M 0 . Fig¬ 
ure [1] compares this model of the LMC gaseous disk to 
the observed HI profile, as found from averages of radial 
profiles in four distinct quadrants of the system. The 
observed HI colu mn was taken from t he Parkcs HI Sur¬ 
vey (HIPASS, see lPutman et al.ll2003l ). This comparison 
shows our model provides an excellent match to the cen¬ 
tral HI column of the observed system, and matches the 
radial falloff in two quadrants of the system (northwest 
and southeast, ignoring the heavily star forming 30 Do- 
radus region). In the southwest quadrant, leading into 
the Magellanic Bridge, the observed column predictably 
rises above the radial falloff. The worst agreement is be¬ 
yond 4 kpc on the leading edge (northeast). This is itself 
evidence of a role for ram pressure stripping, especially 
in light of the uninterrupted presence of stars beyond the 
HI truncation in this quadrant. 

Equating our model’s surface density to the observed 
HI surface density ignores other phases of the ISM. How¬ 
ever, dense molecular clouds behave more akin to the 
stellar component in a ram pressure stripping scenario, 
and thus could be thought of as a small error in the stel¬ 
lar disk mass. And while we do not expect ionized HII to 
be a large contribution to the gas profile by mass, we do 
account for HII in our simulations, crudely, v ia a densi ty 
cut a t p = .03 cm' 3 (see, for example, [Rahmati et al.1 
l2fH3h . 


This model also produces a total gas mass and HI di¬ 
ameter in line with th e tight relation found for a survey 
of late-type dwarfs bv iSwaters et al.l (I2002D . The faster 
falloff along the leading edge today is not in line with 
this relation, which further suggests it is not a long-lived 
remnant of a pre-MW LMC but rather induced by RPS, 
RP compression of contours or RP-induced star forma¬ 
tion and feedback. Our model is on the lower mass end 
of all models that obey both the Swaters relation and the 
mass constraint of K98. Thus in section H~6l we consider 
a more gas-rich initial disk profile as well. 

Our LMC disk sits within a large DM halo, which we 
treat as a spherical potential whose center is matched 
with those of the stellar and gaseous distributions. In¬ 
ferring the circular velocity and dynamical mass of the 
LMC is not a straightforward process. Measurements 
of the rotation curve from different stellar and gaseous 
populations yield disparate results, convoluted further by 
imprecise k now le dge of the system ’s orientation and el- 
lioticit v (Ivan der Marel et al.l 120091 section 2.5 for a re¬ 
view). lOlsen et al.l ( 201ll ) found the HI rotation curve 
peaks at 87 ± 5 km/s beyond 2.4 ± .1 kpc. More 
recently, (Ivan der Marel fc Kallivavalill 12014 hereafter 
vdM14) took 3D stellar proper motions and consistently 
fit the bulk motion, orientation and rotation curve of the 
LMC disk, assuming a flat, circularized stellar disk, to 
find a flat rotation curve of91.7±18.8 km/s, correspond- 
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Figure 1. Comparison of the observed HI column density from the Parkes HI Survey (HIPASS, see lPutman et aDbOO /f) to our 
model’s gaseous surface density (red line). Here we plot surface density versus radial distance from the LMC kinematic center 
in the plane of the LMC disk. The inset shows the observed HI surface density map from which our profiles are derived. To 
map the observed two-dimensional line-of-sight (LOS) column’s position to a radius in the disk plane, we assumed the gas’s 
height above the plane was zero and made use of the coordinate systems and transformations described in Section 14.21 Next 
we sampled the observed profile along a dozen rays, equally spaced in angle from 30° to 60° in each quadrant (grey lines with 
dots, white lines on the inset) and then averaged these radial samplings into an observed radial profile (thick black line). The 
central surface density of our model matches the observed profile well, and the radial falloff is in agreement in the northwest and 
southeast quadrants. The leading edge (northeast) shows a stronger radial falloff beyond 4 kpc, itself evidence of ram pressure 
stripping, as the stellar component does not show any such curtailment here. Figure [6] shows the observed stellar profile. 


ing to a dynamical mass of M( 8.7 kpc) = 1.7 x 10 10 Mq. 
Our model produces a rotation curve of 80 km/s, which is 
consistent with vdM14, though admittedly on the lower 
end of their error space. Our dynamical mass is thus 
« 1.4 x 10 10 M q (includes gas, stars and DM). We also 
perform a few simulations with peaks of 90 and 100 km/s, 
though we find no significant effect from this change. We 
emphasize that an accurate rotation speed is of secondary 
importance to our results, since where RPS truncates the 
disk gas does not depend strongly on this value (see Sec- 
tion[3]). For the D M, fo llowing T09, we employ the static, 
spherical model of lBurkertl (119951 ) . 

Pdm(t) = PqTq [(r + ro)(r 2 + 7q)] _1 (3) 

with characteristic density and radius po = 3.4 x 10“ 24 
g/cm 3 and r = 3 kpc, values chosen to produce a rota¬ 
tion curve peaked at 80 km/s and also to closely follow 
the equivalent NFW profile with concentration param¬ 
eter c = 10. Figure [2] shows this rotation curve also 
broken down into mass compone nts, wi th co mparis ons 
to the constraints of K98, vdM14. [Olsen et al.1 (|2011f ). 

In our simulations, as a practical consideration, we im¬ 
plemented a hydrostatic LMC gas halo (LGH) profile, 
whose pressure declines towards the box edge. This setup 
mitigated problems with our initial headwind setup. Our 


inclusion of this LMC gas halo is otherwise immaterial 
since the headwind washes this small gas component 
away well before the LMC’s gaseous disk gets stripped. 
The low density, hydrostatic gas distribution obeys 

T]lgh = GM DM {R)prn p ^ 

3 ki)T 

ulgh = riLGHO ( 7F -) ( - ) (5) 

\J-LGH 0 / \ r LGH0 / 

where tlgho = 2 kpc and tilgho = 1-8 x 10 -3 cm -3 
were chosen to provide ample pressure support near the 
LMC disk but a very low pressure, low density ambi¬ 
ent medium near where the simulation wind first prop¬ 
agates into the domain. T LG ho = Tlgh(tlgho), and 
these expressions lead to a total LMC gaseous halo mass 
of roughly 5 x 1O 6 M 0 , or roughly 1% the LMC gaseous 
disk mass. An appreciably higher fraction of baryons 
could reside in this halo without affecting our results, 
since the material, spread over the halo’s large volume, 
would still be very diffuse, and thus easy to unbind long 
before pericentric passage. 

2 .2. Model for the MW Diffuse CGM 

For the extended, diffuse phases of the circumgalactic 
medium, we adopt a /3-model in density dMakino et al.l 
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To completely specify the thermodynamic state of our 
CGM model we’ll need not only this /3-profile density and 
our universal ideal gas law equation of state, but also a 
temperature profile. The choice of temperature profile 
is of small significance, since simulations both with and 
without gas cooling have found the stripping dynamics 
strongly insensitive t o the headw i nd’s tempera ture and 
pressure (e.g. T09. iRoediger fc Henslerl [200511 . Again 
appealing to cluster observations^ we choose to mo del 
the MW halo temperature with dMakino et al.lfl99 8) 


k b T(r) = 7 


G/im p M(r) 
3 r 


( 8 ) 


where 7 ss 1.5 is found in numerical simulations. For 
the mass of the MW , we follow the classic NFW profile 
(jNavarro et al 1 ITTTotI) with M 2 oo = 10 12 Mq and concen¬ 
tration parameter c = 12 . 


Figure 2. Our model LMC’s rotation curve, broken down 
by mass components. The net rotation curve peaks at 80 
km/s, on the lower end of current observations. Observation- 
ally inferred circular velocities are over-plotted with errors: 
vdM14 estimated the ampl itude of t he circular velocity from 
stellar kinematic s, whereas I Olsen et al.l (120111 1 did so via HI. 
IKim et all (119981 1 also used the HI velocity field along with 
the luminosity of the baryonic components to find an HI mass 
within r = 4 kpc. We emphasize that correctly modeling the 
stellar and gaseous surface densities at the observed HI trun¬ 
cation radius of fts 6 kpc is the chief concern for our RPS 
study of the LMC disk, whereas the rotation curve profile in 
interior regions is of limited importance. 


119981) . motivated by observations of the hydrostatic gas 
filling large X-ray clusters, 


n(r ) = n 0 



( 6 ) 


where no and r c describe a “core” density and radius 
within which the profile is roughly flat, and /3 is a power 
law decay constant app licable at larger radii r > r c . 
iMiller fe Bregma^ (12013H . hereafter MB 13, fit a /3-profile 
to the distribution of ~ 10 6 K diffuse halo gas to match 
OVI and OVII absorption line strengths for 29 X-ray 
sources (primarily AGN). They found best fit parame¬ 
ters for the MW, 


tto = 0 .461°35 4 cm 3 
r c = 0 . 35 t° 2 7 9 kpc 
/3 = 0.7ll°i4 3 


which suggest the MW core region ends at radii well 
below our region of interest (r ~ 50 kpc). Thus for our 
LMC this profile can be approximated as 


n(r) 


n 0 r 


3/3 

C 


7*3/3 


(7) 


to good precision — i.e. our model is more or less a 
simple power law and thus a function of two important 
variables: an amplitude nor 3 ^, and a falloff /3. Extrapo¬ 
lated mass estimates for the CGM will depend most sen¬ 
sitively on /3. Section 14.41 details how we chose various 
parameters in our suite of simulations. 


2.3. Model for LMC Orbit 

Our work relies on a model of the LMC’s orbital his¬ 
tory to produce the CGM headwind’s time-evolving den¬ 
sity and 3D velocity. Such an orbit was generated via 
a backwards-orbit integration for the LMC/MW system 
using methods described in B07. These calculations treat 
the MW as an NFW profile and the LMC as a Plum¬ 
mer profile and ignore the LMC’s major companion, the 
Small Magellanic Cloud (SMC). But such simplifications 
will not have a huge impact on the accuracy for the 
most critical past 100 Myr surrounding perice ntric pas¬ 
sage (|Kallivavalil et al.ll2013HBesla et alJl2007D . 

The initial conditions for the LMC’ s motion are the 
latest proper motion measurements of iKallivavalil et al.1 
(|2013f l (hereafter K13). For our fiducial orbital sce¬ 
nario, we adopt an orbit with the pericentric velocity 340 
km/s at a distance of 48.1 kpc from the Galactic Cen¬ 
ter. This fiducial orbit was calculated assuming a MW 
with M v ; r = 1 x 1O 12 M 0 , c v ir = 9.86 and R v u = 261 kpc 
with MMw.disk = 6.5 x 1O 1 O M 0 , as in K13. Our fidu¬ 
cial model involves an LMC Plummer sphere with mass 
1.8 x 10 u Mq with dynamical friction calculated also as 
in K13. This dynamical friction implies that uncertainy 
in the MW/LMC mass ratio results in uncertainty in the 
LMC’s orbital history, which we explore in Section [5] 

As described in the introduction to this paper, the 
LMC’s orbital period is an appreciable fraction of a Hub¬ 
ble time, implying either a first-infall scenario or an orbit 
in which the last pericentric passage occurred much far¬ 
ther from the Galactic Center. Thus its recent pericentric 
passage with a speed of « 340 km/s involved a ram pres¬ 
sure headwind of unprecedented strength for the LMC. 
A plot of this evolving ram pressure appears in Figure 01 

2.4. Defining the Truncation Radius of the LMC’s 
Leading Edge 

To compute the observed leading-edge truncation ra¬ 
dius, we took the HIPASS HI column density map and 
re-projected the HI column into our LMC line-of-sight 
(LOS) point-of-view (POV) (see Sections 14.21 and 14.51) . 
We then extended a dozen rays radially from the kine¬ 
matic center of the LMC out to the leading edge, which 
we define as the second quadrant of a cartesian frame 
centered on the LMC kinematic center where the y -axis 
is tangent to the line of constant declination passing 








































6 


Salem et al. 


through the origin (see Section ITT]) . Table 2 summarizes 
position and orientation values relevant to the present 
work. We then find the radial distance at which the col¬ 
umn drops below 10 19 cm -2 for the final time (Section 
I4.5l iustifies this choice). Finally, we fit an ellipse to these 
values whose axis-ratio is constrained to correspond to a 
circle in the LMC disk plane, minimizing the sum of the 
square of the residuals (SSR). The best fit ellipse has a 
corresponding radius of 6.2 kpc in the LMC disk plane. 
The uncertainty in this value greatly exceeds any number 
implied by the goodness of fit of our constrained ellipse or 
concerns over the LMC disk plane orientation. Rather, 
this uncertainty is set by the floccnlent nature of the 
observed HI distribution, with super bubbles and large 
whorls induced by star formation, feedback and the ram 
pressure dynamics itself. By eye, we judge a conservative 
(large) standard deviation of this value to be ~ .25 kpc. 

3. INFERRED MW CGM DENSITY FROM ANALYTIC 
ARGUMENTS 

In this section we take the model outlined in Section 
[2] and use analytic methods to demonstrate how RPS 
considerations coupled with the LMC’s observed leading 
edge truncation radius, Rt , can lead to tight constraints 
on the MW CGM gas density at r ~ 50 kpc. Section 
[5] embarks on a systematic enumeration of our model 
uncertainties, which we ignore in the present section for 
the sake of clarity. 

G72 explored how the ram pressure experienced by 
matter plunging into a denser cluster environment can 
alter the extent of a galaxy’s gas disk by equating the 
ram pressure headwind experienced by the galaxy with 
the restoring pressure of the disk, the latter being propor¬ 
tional to both the gaseous and stellar surface densities, 

pv 2 = 2ttG(E*(R) + a£ g )£ g (I?) , (9) 

where p is the physical density of the ambient medium 
the galaxy passes through, in our case, the MW’s dif¬ 
fuse, gaseous halo; v is the speed of the galaxy’s motion 
relative to this gas; E*(i?) is the stellar surface density 
at cylindrical radius R from the galaxy’s center; E g {R) 
is the surface density of disk gas, and a £ [0,1] repre¬ 
sents the degree to which the gas experiences its own self¬ 
gravity during stripping. For MW-sized systems falling 
into cluster environments - the scenario explored in a ma¬ 
jority of the RPS literature - a = 0 has been employed, 
since the gas fraction is quite low in these systems. How¬ 
ever, for gas-rich late-type dwarves like the LMC, this 
restoring force cannot be ignored. We choose a = 1.0, 
i.e. full se lf-gr avity, a choice well-supported by analysis 
in Section ITHI 

According to this model, (hereafter the GG model), re¬ 
gions of the disk where the ram pressure exceeds the grav¬ 
itational restoring force would experience a rapid loss of 
gas. For a simple exponential disk profile, this restor¬ 
ing pressure falls off monotonically, implying there ex¬ 
ists a truncation radius beyond which gas is completely 
removed. This simple model thus provides a map be¬ 
tween the wind pressure and a truncation radius. For real 
systems where there are considerations of the stripping 
timescale, nearby extraplanar material, and the effects 
of viewing angle, the truncation radius predicted by GG, 
Rgg, may not match the observed truncation radius, Rt, 
precisely. Thus for our simulation work, we will consider 


an extended model 

logio(n) = log 10 [n G G(Rt)] + D (10) 

where 

27tG E gas (I?t)[E*(I? t ) + E gas (I?t)] ^ 

nGG{Rt) = - n - ■ ( 11 ) 

prrip v 

Here n is the actual CGM halo density, whereas uqg is 
the prediction of the GG model. D is an offset term 
meant to capture any observed discrepancies between 
the two. In the present section, where we explore the 
model’s capabilities using simple analytic calculations, 
we set D = 0. In section [5] we use data from a large 
suite of 3D hydrodynamic simulations to better describe 
the offset term and its effect on the implied MW CGM 
density. 

The LMC’s velocity during its most recent pericentric 
passage is well constrained (K13), as is its gaseous and 
stellar surface density profiles. Within the GG model, 
this leaves only the ambient CGM’s density, p, as a free 
parameter, which given the rapid stripping timescales of 
R05 and T09, we may interpret as the density at peri- 
galacticon. Thus GG together with the LMC’s observed 
HI and stellar extent imply a localized measurement of 
MW CGM density at LMC perigalacticon. 

To explore this approach, we sampled the /3-profile 
(Equation [5]), selected a fiducial LMC orbital history 
(defined in Section 12.31) . used the implied halo density 
at perigalacticon and orbital speed to find the peak ram 
pressure, and then applied the GG model of Equation [G] 
to find a predicted truncation radius for the LMC disk’s 
leading edge. We then explored which CGM models 
matched the observed HI truncation radius. 

5,000 /3-profiles were randomly selected, drawing the 
three parameters from independent, uniform distribu¬ 
tions: the core density no £ [0,1] cm -3 , r c £ [0,1.2] 
kpc and /? £ [0,1.4]. The choice of these bounds were 
guided by MB13 (see Figure [3]). We selected the fiducial 
orbit of K13, which at perigalacticon had r p = 48.1 kpc 
and v p = 340. km/s. In section [5] we explore how the 
spread in allowed orbits affects our results. With this 
speed and distance from the Galactic Center, Equation 
[G] provides a mapping to the density at perigalacticon, 
rip, and peak ram pressure, P p , for each halo model. We 
then apply GG via Equation [9] solving numerically for 
the gaseous disk’s truncation radius, Rgg, assuming the 
LMC gas and stellar distribution follow the model out¬ 
lined in Section m via Equationl and 2. Thus for each 
halo model, there exists an implied truncation radius, 
Rgg, which can be compared to the observed value. 

Figure [3] shows the distribution of our 5,000 halo mod¬ 
els in parameter space, both for the full /3-profile, in 
terms of n o, r c and /3 and for the approximation of Equa¬ 
tion [3 valid at LMC distances and described by /3 and 
nor 3/3 . We highlight models that successfully fall within 
1-, 2- and 3-er of the observed truncation radius. The top 
four panels of Figure [3] show that this constraint places 
a strong upper limit on /3, for reasonable core density 
choices (governed by no and r c ), and for a fixed core pro¬ 
file, selects a narrow range of allowed /? exponents. This 
reflects how much leverage our measurement at ~ 50 kpc 
from the Galactic Center has, which strongly constrains 
any simple exponential halo profile’s falloff. In contrast, 
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Figure 3. Parameter space for a /3-profile model of the MW’s diffuse CGM. The core density, no, core radius, r c and exponential 
falloff, /?, completely specify the profile, described by Equation [6] The scatter plots show 5,000 model choices, randomly 
selected from independent, uniform distributions of the three parameters. The GG model (Equation [9]) was then used along 
with our LMC model’s stellar and gaseous distributions to determine a predicted RPS truncation radius of the LMC’s leading 
edge. Models where this value is within l-er of the observed radius (6.2 ± .25 kpc) are shown in magenta, models within 2-<r 
in cyan and 3-cr in orange. These cuts separate cleanly only in the bottom-right panel, where we’ve plotted the “parameters” 
/3 vs. nor®' 9 , motivated by the simplified Equation^ valid for r/r c /§> 1. Over-plotted on each panel, we show the 1-, 2- and 
3-a confidence intervals for the constraints found by MB13 (their Figure 4), from X-Ray emission and quasar absorption lines. 
Towards the end of our investigations, MB15 published tighter constraints on this parameter space, shown in the bottom 
contour plot, with our analytic results again over-plotted. Figure [14] reproduces this bottom panel, but with the results of our 
3D hydrodynamic simulations and full error analysis folded in. 
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Quantity 

Value 

Unit 

Source 

Description 

a?o 

78.76 

deg 

vdM14, Tab. 1, col. 3 

CM position in celestial 


-69.19 

deg 

vdM14, Tab. 1, col. 3 

coords 

D 0 

50.1 

kpc 

Freedman et al. f20014 


v x 

-453.8 

km/s 

K13 

Orbital velocity in the 

Vy 

54.4 

km/s 

K13 

LOS frame (Sun at rest) 

V z 

-262.2 

km/s 

K13 


e 

139 

deg 

vdM02 

Disk pos/inc angles and 

i 

34.7 

deg 

vdM02 

angular change rate 

di/dt 

0 

deg/Gyr 

vdM14 


^sys 

262.20 

km/s 

(calculated) 

Recession velocity 

Vt 

457.05 

km/s 


Transverse speed 

0t 

83.16 

deg 


Position angle of Vt 


Table 2: Summary of LMC COM position and velocity as well as the disk plane’s orientation and solid body motion. The top 
three quantities describe the LMC COM’s position, with RA and DEC (ao and /3o) and the distance from the Sun to the LMC, 
Do. The next three quantities describe the COM’s proper motion in the LOS frame (See Section f4. 21 . The third set describe 
the LMC disk plane’s orientation on the sky via the position angle of the line of nodes O, the plane’s inclination angle with 
respect to the sky plane i, and the rate of change of this inclination angle di /ctt which alters the appearance of internal motion 
within the LMC (see Section [4.51) . The final set of quantities uses Equation I A2I to translate the COM into a form amicable to 
the velocity corrections applied in Section T4. 5 1 
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the LMC’s orbit does not afford us a strong probe of 
the innermost region of the CGM, as evidenced by the 
nearly uniform spread in no - r c space. The /3 - nor^f 
space of the approximation in Equation [7] provides the 
cleanest separation, and selects out the narrowest range 
of /3 for a given core density. Also over-plotted in the 
top four panels of Figure [3] are the 1-, 2- and 3-cr con¬ 
tours from MB13’s work analyzing direct observations of 
~ 10 6 K gas. These results suggest our approach pre¬ 
dicts consistently lower /3-values (i.e. higher density at 
large radii) than their results, though their is good agree¬ 
ment for the smallest r c and no core parameters shown. 
Towar ds the end of our investigations. lMiller &; Bregmanl 
(j2015H (hereafter MB15) published tighter constraints on 
the allowed halo parameters, narrowing the results to /3- 
profiles with small cores ( nor 3/3 ~ .01). The final panel of 
Figure [3] plots the agreement of our analytic RPS model 
with the observed truncation radius with MB15’s up¬ 
dated error space superimposed. We again find rough 
agreement, though the 1-sigma contours just barely fail 
to overlap. Averaging the density at LMC pericenter 
from these samples, we find n(48 kpc) « 8 x 10 -5 cm -3 . 
Note our results in these plots represent only an analytic 
toy model. We reproduce this panel in Figure HU with 
the results of our hydrodynamical simulations and full 
error analysis. 

This crude approach applies best to a smooth, az- 
imuthally symmetric LMC disk of low gas-fraction, fac¬ 
ing a strong, constant, head-on wind (aligned with the 
disk’s angular momentum axis). In reality, a wind of vari¬ 
able strength arrives at a varying angle to an LMC disk 
riddled with large scale inhomogeneities and a significant 
velocity dispersion. On a localized level, the GG model 
itself comes into question, as the details of momentum 
transfer between the wind and gas involve shock fronts 
propagating in front of and behind the contact point be¬ 
tween the wind and disk, making this as much a thermal 
process as a mechanical one. In addition, the flocculent, 
multiphase ISM is threaded with magnetic fields and in¬ 
fused with photons and energetic particles, complicating 
notions of its tensile strength and role as a monolithic 
impediment to the wind’s motion. Since pericentric pas¬ 
sage happens quickly, and has occurred less than 100 Myr 
ago, the validity of using the ram pressure at this epoch 
depends on the timescale of RPS. The GG argument’s 
only (implicit) timescale is the crossing time of the wind 
through the galaxy, set by the windspeed and the disk’s 
scale height, which may play no role in setting the actual 
stripping timescale. 

Despite this, many numerical simulations have con¬ 
firmed GG’s applicability i n a host of stripp ing scenar¬ 
ios. In particular. iRocdigcr fc Henslerl ( 2005 1. hereafter 
R05, simulated ram pressure in two dimensions, and 
looked across many wind strengths and vertical disk scale 
heights, and found the timescale of fast stripping to oc¬ 
cur ~ 20 — 200 Myr, followed by longer, slower phases 
of stripping, the final of which was driven by viscous 
processes. The results were only strongly dependent on 
the wind properties (not the disk scale height), most no¬ 
tably its ram pres sure, and i n good agr e emen t with GG. 
In their followup. iRoediger fc Briiggenl (|2006H . hereafter 
R06, performed fully three-dimensional runs, showing 
that disk orientation does not strongly alter the stripping 
extent except for extreme angles (beyond ~ 60°). Be¬ 


yond this, T09 explored stripping for a host of ram pres¬ 
sures both with and without radiative cooling, and found 
agreement to within 10% of GG. These explorations sug¬ 
gest the model of Equation [9] is a powerful, predictive 
tool relevant to our studies. In particular, the LMC’s re¬ 
cent orbital history involves a disk galaxy having passed 
through a period of strong, peak ram pressure over the 
last ~ 100 Myr, with an orientation angle between wind 
and disk vector well within the ~ 60° limit of R06. 

4. INFERRED MW CGM DENSITY FROM SIMULATIONS 

Motivated by the analytic proof-of-concept in Section 
[U we now present three-dimensional, global galaxy simu¬ 
lations of an LMC-analog isolated disk galaxy with stel¬ 
lar, gaseous and dark matter components, subjected to a 
time-varying headwind meant to mimic the diffuse CGM 
of the MW as the dwarf swoops through its recent peri¬ 
centric passage. The chief goal of this work is to corrob¬ 
orate and extend the model above, though other useful 
information will arise from these investigations, includ¬ 
ing a look at how RPS influences the LMC’s HI velocity 
structure, and an upper bound on how much HI RPS 
can redistribute from the LMC disk to the MS. Sections 
14.11 through 14.41 describe the numerical implementation, 
wind setup and halo model selection. Some of this ma¬ 
terial is presented mostly for the sake of reproducibility, 
and the casual reader may wish to skip details, particu¬ 
larly Section 1431 Section l4~5l provides the results of these 
simulations, building towards a thorough understanding 
between the CGM density and the LMC’s leading edge 
truncation radius. 

4.1. Numerical Implementation 

For the present work, we employe d the Eulerian hydro¬ 
dynamics code Enzo, described in iBrvan et all ( 2014 1. 
Our runs are fully three-dimensional and make use of 
Enzo ’s robust ZEUS hydro method dStone fc Normanl 
I1992T ). Self gravity of the gas is also included. 

Our LMC sits within a 60 kpc-wide cubic box, with a 
256 3 base grid resolution. One of Enzo’s main strengths 
is adaptive mesh refinement (AMR), which uniformly re¬ 
solves the entire simulation region on a course grid but 
provides higher resolution, “refined” sub-grids as needed 
in regions where the dynamics grow complex. In our 
work, gas density a factor of four above the background 
density triggers higher mesh refinement, on up to three 
additional levels, refining by a factor of two each time. 
This leads to an effective resolution of 30 pc throughout 
the entire disk, or less than a tenth of the disk’s vertical 
scale height. 

We use a “wind tunnel” setup for our simulations, 
working in a frame where the LMC sits at rest within 
our simulation domain, and the MW Halo thus becomes 
a “headwind”, produced from the boundary of the sim¬ 
ulation box. This setup reduces computational costs in 
a few ways: first it guarantees the densest, most-refined 
disk material is at rest in our simulation frame, which 
prevents relative motion of this gas from limiting the 
simulation’s time stepping; and second it alleviates the 
need to simulate the entire volume of the Milky Way halo 
through which the LMC moves. We place our galaxy 20 
kpc in each dimension from the corner of our box desig¬ 
nated as the origin, where the halo wind first propagates 
inward. Our initial distribution of baryons follows the 
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idealized LMC model outlined in Section l2Jl and is con¬ 
sistent with a three-dimensional, time-dependent exten¬ 
sion of the analytic analysis of Section [3j Our gaseous 
disk is evolved in isolation, before any ram pressure head¬ 
wind arrives, for 1 Gyr, to allow transient oscillations 
in the structure to die away. We then launch the wind 
to mimic the LMC’s orbit, beginning 1 Gyr before the 
present day, long before the ram pressure’s impact grows 
significant. 

Our simulations include self-gravitating gas, represent¬ 
ing the atomic and ionized components of the LMC’s 
ISM and the MW gaseous halo. We elect to shut off ra¬ 
diative cooling and star formation and feedback in this 
simulation. This choice was motivated by an acknowl¬ 
edgement that the current generation of feedback mod¬ 
els fail to produce an accurate ly pressurized ISM (al¬ 
though see iHookins et all 12014b , and thus inclusion of 
this physics may lead to an ISM with erroneously low 
tensile strength, thereby overestimating its susceptibility 
to stripping. Stars and dark matter are included in a 
static sense: we compute the gravitational acceleration 
caused by the stellar disk and LMC DM halo distribu¬ 
tions discussed in Section IQ and include these time- 
invariant, spatially varying acceleration vectors into the 
simulation. The potential of the MW and forces associ¬ 
ated with our accelerating, wind-tunnel reference frame 
are ignored, as our system’s dynamics of interest lie well 
within the likely LMC tidal radius, which vdM14 found 
to be r = 22.3 ± 5.2 kpc. To cleanly separate the role 
of RPS from other LMC-specific dynamics, we do not 
include the SMC, which has a strong gravitational influ¬ 
ence on the LMC. 

4.2. Coordinate Systems 

Foll owing the work of Ivan dcr Marel fe Cionil (|200l! ) 
and Ivan der Marel et al.l (120021) (hereafter vdM02), we 
define the following four important Cartesian reference 
frames used in the present work. 

Line of Sight (LOS) Frame: — A 3-D cartesian coordinate 
system whose origin is the “center” of the LMC. For our 
simplistic model of the LMC, this corresponds to the 
center of its circular gaseous disk’s initial distribution; its 
static, circular stellar potential; and its static, spherically 
symmetric dark matter potential. The i-axis lies anti¬ 
parallel to the right ascension (RA, a, at the origin of 
the system), the y-axis parallel to the declination (d) and 
the f-axis pointed towards the solar neighborhood (i.e. 
the observer). From this frame we can easily collapse 
our 3-D data into 2-D projections of the LMC that can 
be overlaid on and compared to real observations with 
relative accuracy. 

LMC Frame: — A n other frame first devised in 
Ivan der Marel &; Cionil (120011) centered on the LMC, 
with the vertical z'-axis aligned antiparallel to the 
galaxy’s angular momentum vector, L. In the present 
model, L is perfectly aligned for the stars, gas and DM. 
This new frame is defined in terms of two rotations out of 
the LOS frame: first a counterclockwise rotation about 
the z (line-of-sight) axis by an angle 9, and second a 
clockwise rotation about the new i'-axis by an angle i. 
Here i corresponds to the inclination angle of the LMC 
disk out of the LOS plane: the plane on the sky passing 


through the LMC center whose normal vector is aligned 
with the line-of-sight. And 9 corresponds to the angle 
between the horizontal LOS axis and the line of nodes, 
which represents the intersection of the LMC disk plane 
and the LOS x-y plane. For these angles, we use the 
results of vdM02 obtained from analysis of carbon star 
kinematics, setting i = 34.7° and O = 9 — 90 = 139.9°. 

Simulation (SIM) Frame: — Working in the LMC frame 
is an attractive option, since height above and below the 
disk plane is aligned with the z'-axis and all of the disk 
gas’ circular motion is in the x'-y' plane. This not only 
simplifies how enzo generates the initial conditions but 
also aids in the solver’s accuracy, since angular momen¬ 
tum is better preserved when rotational motion is aligned 
with the grid axis. However, one downside of the frame 
as defined by Ivan der Marel fc Cionil ( 200D is that the 
headwind experienced by the LMC would have veloc¬ 
ity components that switch sign over the course of the 
LMC’s orbit. This means a wind being blown in from 
one boundary face of our simulation box would eventu¬ 
ally switch to being blown in from the opposite bound¬ 
ary face — a complication that presents a formidable 
software challenge. We can easily circumvent this issue 
by rotating the LMC frame 100° clockwise about the z'- 
axis, retaining the LMC frame’s attractive features. All 
our present simulations were run in this new SIM frame. 

Galactic Frame: — This last frame of interest is the 
only one not centered on the LMC but rather on the 
solar nei ghborhood, and where the LMC orbit data from 
iBesla et al.l (2012) is most naturally defined. Here the 
z MW -axis is aligned with the MW’s angular momentum 
vector, and the f MW axis points away from Sgr A*. Here 
the following rotation matrix pulls a vector in LOS co¬ 
ordinates into this frame 


MW ' 


'0.11638 -0.98270 -0.14410' 

MW 

= 

0.57132 -0.05244 0.81905 

MW 


-0.81243 -0.17765 0.55533 


( 12 ) 

For the sake of both clarity land reproducibility, Table 
3 displays important vector quantities used in this work, 
instantiated in our four coordinate frames. 

4.3. Simulation Boundary Wind 

Our simulated LMC stands still in our simulation box, 
with its angular momentum axis aligned with the carte¬ 
sian grid. To capture the hydrodynamical effects of its or¬ 
bit through the MW, we launch a wind into the box from 
three orthogonal, adjacent sides of the domain. The op¬ 
posite three sides are set to “outflow” conditions, to allow 
the downstream gas to escape the domain. The density, 
temperature and three-dimensional velocity components 
of this wind evolve over time to mimic the headwind ex¬ 
perienced by an orbiting LMC in its final Gyr before the 
present. 

We take our fiducial LMC orbit from K13 and assume 
the orientation of the LMC angular momentum axis is 
constant in a Galactic frame. Neither this particular or¬ 
bital scenario nor this constant orientation assumption 
will strongly influence our results (see Section [5]). This 
orbital data supplies us with the time-evolving three di- 
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Quant 

LMC 

SIM 

LOS 

Galactic 

£ lmc 

[1., O., 0.] 

[-0.174, -0.985, 0.0] 

[-0.656, -0.755, 0.] 

[0.665, -0.335, 0.667] 

y LMC 

[ 0 .. 1 ., 0.] 

[0.985, -0.174, 0.0] 

[0.620, -0.539, -0.569] 

[0.684, -0.083, -0.724] 

*LMC 

[0., 0 ., 1.] 

[0., 0 ., 1.] 

[0.430, -0.373, 0.822] 

[0.299, 0.939, 0.174] 

*LOS (WEST) 

[-0.656, 0.620, 0.430] 

[0.725, 0.538, 0.430] 

[1.. 0 ., 0.] 

[0.116, 0.571, -0.812] 

3/LOS (NORTH) 

[-0.755, -0.539, -0.373] 

[-0.400, 0.837, -0.373] 

[ o., 1 ., 0.] 

[-0.983, -0.052, -0.177] 

*LOS (-AWAY) 

[0., -0.569, 0.822] 

[-0.561, 0.099, 0.822] 

[0., 0 .. 1.] 

[-0.144, 0.819, 0.555] 

Angle between disk 
and cc-y-plane 

0.0° 

0.0° 

34.7° 

76.3° 

L 

[0., 0 ., -1.] 

[0., 0 ., -1.] 

[-0.430, 0.373, -0.822] 

[-0.299, -0.938, -0.174] 

v Sun 

[-72.0, -18.6, 240.3] 

[-5.85, 74.1, 240.3] 

[138.9, -25.4, 208.2] 

[11.1, 251.2, 7.25]* 

V LMC ,unc 

[256.7, -161.7, -430.8] 

[-203.8, -224.7, -430.8] 

[-453.8, 54.4, -262.2] 

[-68.5, -476.8, 213.4]** 

V LMC = 

V LMC,unc "I" v Sun 

[184.6, -180.3, -190.5] 

[-209.6, -150.5, -190.5] 

[-314.8, 29.03, -53.951] 

[-57.4, -225.6, 220.7] 

Angle: v LMC , L 

53.571° 

53.571° 

53.571° 

53.571° 

Angle: v wind , L 

126.429° 

126.429° 

126.429° 

126.429° 


* See discussion in Section 4.1 of Ivan der Marel et al.l 120121) . ** Value from K13 

Table 3: Various vector quantities expressed in our four cartesian coordinate frames defined in Section roi Bold 
quantities represent observed inputs or quantities expressed in their most natural coordinate frame. Here L denotes the LMC 
disk’s angular momentum unit vector. vs un is the Sun’s velocity with respect to the Galactic Center (GC). vlmc.uiic is motion 
of the LMC COM with respect to the Sun, whereas the “corrected” vlmc is relative to the GC. This last quantity is the most 
useful in constructing our ram pressure headwind, which is denoted here at present day as v w i n d = — vlmc- 


mensional headwind velocity and the evolving radial dis¬ 
tance of the LMC from the MW center which we can 
transform into our simulation’s frame of reference with 
the aid of Section 14.21 Equations [6] and [3] map the ra¬ 
dial distance to the wind’s density and temperature, thus 
completing its description. 

Our simulated wind is an orchestrated attempt to make 
fluid inflow from three orthogonal boundary surfaces 
of our domain to mimic gas in hydrostatic equilibrium 
within the MW’s gravitational potential, as experienced 
by a dwarf system whose local rest frame is non-inertial. 
Thus achieving a good match between model and real¬ 
ity throughout our simulation box is a non-trivial mat¬ 
ter that needs to be verified “experimentally”. Figure 0] 
shows how this wind evolves over time for our “fiducial” 
simulation M-1.2 (see Section HU). The figure displays 
not only the input wind given to our simulation software 
but also a measurement of the wind probed directly in 
our simulation box halfway between the box’s origin and 
the LMC kinematic center, upstream of the galaxy. This 
probe demonstrates excellent fidelity between the theo¬ 
retical wind evolution and the wind actually experienced 
by the simulated disk. The initially wide disparity be¬ 
tween modeled and simulated velocity components sig¬ 
nifies the presence of an initial “adjustment” shock, and 
does not affect our results since the ram pressure at this 
early stage is far too weak to alter the disk’s gas struc¬ 
ture. The fidelity of the propagating wind was also veri¬ 
fied for larger distances from the box edge in a simulation 
devoid of an LMC itself. Figure 0] also illustrates how the 
ram pressure is far stronger at the LMC’s recent pericen¬ 


tric passage than at earlier times in its orbital evolution, 
suggesting its current HI extent is a localized probe of 
the hydrodynamics at r p ~ 48 kpc. 

The changing components of the wind’s three dimen¬ 
sional velocity in Figure 0] demonstrate how the orienta¬ 
tion between the LM C disk and the headw i nd ch anges 
measurably over time. iRoediaer fc Briiggen ( 200 6) sim¬ 
ulated ram pressure stripping for a variety of wind-disk 
inclination angles, and showed that the gas disk rapidly 
strips to the truncation radius predicted analytically by 
Gunn & Gott in all but the the most edge-on orienta¬ 
tions. We explore the role of wind orientation angle fur¬ 
ther in Section 031 

4.4. Simulation CGM Model Selection 

We seek to constrain the CGM density at radii cor¬ 
responding to the LMC’s recent pericentric passage at 
r « 48 kpc using the current observed HI disk extent. 
As discussed in Section 031 ram pressure is strongest at 
the LMC’s most recent pericentric passage, suggesting 
that the resulting HI truncation radius will be mostly 
sensitive to the local halo density at pericenter, n p = 
n(r ss 48 kpc). It is this quantity we systematically vary 
in our simulations. 

However, our simulations must also address what role 
the broader CGM density profile plays for a fixed pericen¬ 
tric density. If the Gunn & Gott picture of ram pressure 
stripping holds precisely, and the process occurs instan¬ 
taneously, then indeed the density at exactly pericenter, 
for a given orbital model, is all that matters. But if the 
stripping time scale, chaotic gas dynamics, oscillations 
within the disk, and the history of the ram pressure faced 
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Time [ Myr ] 

Figure 4. Comparison of intended wind within our simula¬ 
tion box (dashed lines) to the wind actually measured within 
the simulation over time (solid, colored lines). This wind de¬ 
pends on both the LMC orbit and chosen MW gaseous halo 
model. Shown here is our fiducial, M-1.2, simulation. The 
simulation measurement shown here was taken halfway be¬ 
tween the corner of the simulation box where the wind is 
introduced and the center of the LMC. We also tested the 
fidelity of the wind as it propagated further into the simula¬ 
tion domain by running a simulation devoid of an LMC, with 
similarly well-matched results. 

by the disk influence the final outcome, the results may 
become sensitive to the broader picture of gas in the halo. 

While our CGM density model has three parameters, 
no, r c and /3, as explained in Section [2721 our LMC effec¬ 
tively experiences a simpler profile with two parameters: 
an amplitude and an exponential falloff /3. Thus 

while the chief goal of our parameter search will be to 
understand the role of n p , a secondary goal will be to 
explore a few combinations of n$r^P and /3 for a range of 
pericentric densities. 

We chose three representative amplitudes: £ 

{.01, .079, .5} cm -3 kpc 3 , designated as lo-, mid- and hi- 
core simulations, respectively. The mid-core runs were 


given no = .46 cm' 3 and r c = .35 kpc, corresponding to 
the maximal likelihood values of MB13. The other two 
pairs of no and r c were chosen to roughly span two or¬ 
ders of magnitude in the core constant. For each of these 
three scenarios, we then selected /3 chosen to fix n p to 
a wide range of values from ~ 1CP 7 - 10 -2 cm -3 , with 
a denser sampling around the optimal parameter choices 
for matching the observed truncation radius. 

Table 4 summarizes these simulation choices, grouped 
by core constant. The runs follow a standard naming 
convention with a capital letter denoting the core con¬ 
stant level followed by a dash and the halo gas density 
at LMC pericenter in units of 10 -4 cm -3 . For instance, 
the medium-core run with ss .079 and a pericentric 
halo density of 1.2 x 10 -4 cm -3 is denoted M-1.2; the 
high-core density run with pericentric density 7.96 x 10 -5 
is labeled H- . 8; etc. A few miscellaneous runs were also 
thrown in with widely varying halo parameters. Finally, 
we explored the role of the LMC wind orientation by per¬ 
forming a series of runs with varying halo profiles where 
the wind was launched face-on at the galaxy, along its an¬ 
gular momentum axis, denoted with the prefix FL-, FM- 
and FH- for low, medium and high core halos, and some 
runs at a fixed 30-degree inclination angle for the wind, 
denoted IL-, IM- and IH-. Figure [5] also shows each 
simulation choice graphically in parameter space, atop 
the same analytic models of Figure [3] From this over¬ 
lay, it’s clear that our densest sampling of models feature 
systematically lower /3’s (and thus, systematically higher 
pericentric densities) than the optimal values predicted 
by the toy model in Section [3] We discuss and quantify 
this offset in Section [5] 

The simulations described here all make use of our 
standard LMC model, described in Section 12.11 While 
GG predicts that the choice of DM halo has a negligible 
effect on the gas stripping radius, our choice of stellar 
and gaseous disks, in particular their surface densities 
as a function of radius, will strongly impact results. We 
consider uncertainties in the stellar disk via our analytic 
GG model in Section [5] For the gaseous component, we 
ran a smaller suite of simulations with a “gas-rich” LMC 
disk, with the results and implications discussed in Sec¬ 
tion 14.61 

4.5. Simulation Results 

The following analysis made use of yt dTurk et al.l 
l20H for simulation da ta red uction, Astropy 
( Astronv Collaboration et al.l I2013D for observed 
data reduction and map projections, and emcee 
dForeman-Mackev et al.1 12013 ) for MCMC analysis, all 
community-developed, Python packages. 

4.5.1. HI Column Distribution and Inferred Truncation 

Radius 

We begin by describing our fiducial simulation, M-1.2, 
a mid-core halo run with n p = 1.2 x 10 -4 cm' 3 , which 
produced one of the most successful matches between 
the simulated and observed HI truncation radius on the 
disk’s leading edge. The LMC gaseous disk is, at first, 
completely unperturbed by the low density headwind the 
galaxy experiences 1 Gyr ago, at our run’s start. At this 
point, the truncation radius predicted by Gunn & Gott, 
i?GG, lies well beyond the disk’s initial extent. The LMC 
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Lo Core (blue) 

— »„r y .01 




71 0 

r c 

0 

p P 

Total Mass 


r gg 

R t 

Name 

0.139 

0.100 

0.403 

9.04e-14 

1.18e+10 

7.95e-05 

6.09 

6.56 

L-0.8 

0.139 

0.100 

0.396 

1.03e-13 

1.37e+10 

9.06e-05 

5.97 

6.39 

L-0.9 

0.139 

0.100 

0.391 

1.13e-13 

1.52e+10 

9.93e-05 

5.89 

6.17 

L-1.0 

0.139 

0.100 

0.385 

1.26e-13 

1.73e+10 

l.lle-04 

5.78 

6.05 

L-l.l 

0.139 

0.100 

0.381 

1.36e-13 

1.88e+10 

1.20e-04 

5.72 

5.98 

L-1.2 

0.139 

0.100 

0.376 

1.49e-13 

2.09e+10 

1.31e-04 

5.63 

5.81 

L-1.3 

0.139 

0.100 

0.372 

1.61e-13 

2.28e+10 

1.41e-04 

5.57 

5.71 

L-1.4 





Mid Core (green 

) — « .079 



71 0 

r c 

0 

p P 

Total Mass 

np 

r gg 

Rt 

Name 

0.46 

0.35 

1.00 

2.Ole-16 

1.89e+07 

1.77e-07 

12.6 

No Effect 

VHiBeta 

0.46 

0.35 

0.71 

1.46e-14 

1.05e+09 

1.29e-05 

7.82 

6.69 

M-0.13 

0.460 

0.350 

0.680 

2.27e-14 

1.68e+09 

2.00e-05 

7.40 

6.72 

M-0.2 

0.46 

0.35 

0.64 

4.10e-14 

3.22e+09 

3.61e-05 

6.83 

6.83 

M-0.36 

0.460 

0.350 

0.591 

8.45e-14 

7.22e+09 

7.45e-05 

6.15 

6.47 

M-0.75 

0.460 

0.350 

0.586 

9.10e-14 

7.85e+09 

8.02e-05 

6.08 

6.56 

M-0.8 

0.460 

0.350 

0.578 

1.02e-13 

8.97e+09 

9.03e-05 

5.97 

6.56 

M-0.9 

0.460 

0.350 

0.571 

1.14e-13 

1.01e+10 

1.00e-04 

5.88 

6.49 

M-1.0 

0.46 

0.35 

0.57 

1.15e-13 

1.02e+10 

1.02e-04 

5.87 

6.37 

M-1.0-II 

0.460 

0.350 

0.565 

1.24e-13 

l.lle+10 

1.09e-04 

5.80 

6.43 

M-l.l 

0.46 

0.35 

0.559 

1.36e-13 

1.23e + 10 

1.20e-04 

5.72 

6.26 

Fiducial (M-1.2) 

0.46 

0.35 

0.553 

1.48e-13 

1.36e+10 

1.31e-04 

5.64 

6.13 

M-1.3 

0.46 

0.35 

0.548 

1.60e-13 

1.48e+10 

1.41e-04 

5.57 

6.02 

M-1.4 

0.46 

0.35 

0.544 

1.69e-13 

1.59e+10 

1.49e-04 

5.52 

6.18 

M-1.5 

0.46 

0.35 

0.539 

1.82e-13 

1.73e+10 

1.61e-04 

5.45 

5.64 

M-1.6 

0.46 

0.35 

0.535 

1.93e-13 

1.85e+10 

1.70e-04 

5.40 

5.89 

M-1.7 

0.460 

0.350 

0.497 

3.39e-13 

3.52e+10 

2.99e-04 

4.90 

4.89 

M-3.0 

0.46 

0.35 

0.20 

2.74e-ll 

6.34e+12 

2.40e-02 

1.26 

Destroyed 

VLoBeta 





Hi Core (red) 

— n 0 r c ~ -5 




n 0 

T’C 

0 

p P 

Total Mass 


r gg 

Rt 

Name 

0.810 

0.800 

0.751 

9.Ole-14 

6.23e+09 

7.96e-05 

6.09 

6.19 

H-0.8 

0.810 

0.800 

0.732 

1.14e-13 

7.98e+09 

1.00e-04 

5.88 

6.12 

H-1.0 

0.810 

0.800 

0.725 

1.24e-13 

8.75e+09 

1.09e-04 

5.80 

6.30 

H-l.l 

0.810 

0.800 

0.718 

1.35e-13 

9.60e+09 

1.19e-04 

5.72 

6.10 

H-1.2 

0.810 

0.800 

0.711 

1.47e-13 

1.05e+10 

1.30e-04 

5.64 

6.21 

H-1.3 

0.810 

0.800 

0.705 

1.59e-13 

1.14e+10 

1.40e-04 

5.57 

6.19 

H-1.4 

0.810 

0.800 

0.699 

1.71e-13 

1.24e+10 

1.51e-04 

5.51 

6.09 

H-1.5 

0.810 

0.800 

0.689 

1.93e-13 

1.41e-|-10 

1.70e-04 

5.40 

6.04 

H-1.7 





Exploratory (purple) 




71 0 

r c 

0 

p P 

Total Mass 

n P 

r gg 

Rt 

Name 

1.20 

0.64 

0.57 

8.44e-13 

7.50e+10 

7.44e-04 

4.12 

4.10 

MaxWind 

0.80 

0.50 

0.69 

7.12e-14 

5.22e+09 

6.29e-05 

6.31 

6.42 

T-6 

1.50 

1.50 

0.97 

7.03e-14 

5.41e+09 

6.22e-05 

6.32 

6.49 

T-6-X 





Face-On Win 

d (yellow dots) 




n 0 

r c 

0 

p P 

Total Mass 


r gg 

Rt 

Name 

0.139 

0.100 

0.381 

1.36e-13 

1.88e+10 

1.20e-04 

5.72 

6.44 

FL-1.2 

0.460 

0.350 

0.586 

9.10e-14 

7.85e+09 

8.02e-05 

6.08 

7.06 

FM-0.8 

0.460 

0.350 

0.559 

1.36e-13 

1.23e+10 

1.20e-04 

5.72 

6.91 

FM-1.2 

0.460 

0.350 

0.535 

1.93e-13 

1.85e+10 

1.70e-04 

5.40 

6.34 

FM-1.7 

0.810 

0.800 

0.718 

1.35e-13 

9.60e+09 

1.19e-04 

5.72 

6.88 

FH-1.2 




30 c 

J Inclined Wind (Shown in Figure 110i 



71 0 

r c 

0 

p P 

Total Mass 


r gg 

Rt 

Name 

0.139 

0.100 

0.381 

1.36e-13 

1.88e+10 

1.20e-04 

5.72 

6.03 

IL-1.2 

0.460 

0.350 

0.559 

1.36e-13 

1.23e+10 

1.20e-04 

5.72 

6.22 

IM-1.2 

0.810 

0.800 

0.718 

1.35e-13 

9.60e+09 

1.19e-04 

5.72 

6.26 

IH-1.2 


Table 4: Our 44 wind tunnel simulations. The top three groups (L-, M-, and H- runs) each have a different, fixed core halo 
density (determined by no and r c ) with various /3 chosen to explore a range of pericentric densities. The fourth group are runs 
with miscellaneous halo parameter choices. The fifth group is a set of runs where the orbital wind was delivered entirely face 
on to the galaxy. The first three columns, no [ cm” 3 ], r c [ kpc ] and /3, denote the MW halo parameter choices that inform the 
LMC wind density. The next three columns list important implications of these parameter choices: Pram, the peak (pericentric) 
ram pressure experienced by the LMC; the total mass of the MW halo, if we were to extrapolate the /3-profile to large radii 
( M(r > r v ir)); and n p , the MW halo density at the LMC’s pericentric passage. Pgg is the truncation radius of the LMC disk 
inferred from the pericentric ram pressure and the initial disk profile. Rt is the truncation radius for the observable HI inferred 
from the simulations (see Section [4.51 . The final column indicates the name of the run (see Section ^.4l for naming conventions). 
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Figure 5. Our halo model’s parameter space, as in Figurc[3j but now with our 44 three-dimensional, hydrodynamic simulations 
over-plotted as stars. For the majority of runs, the halo’s core parameters, no and r c , were held at three distinct, fixed 
combinations designated “lo core” (blue stars), “mid core” (green stars) and “hi core” (red stars). For each of these core 
choices, /3, the exponential falloff, was incrementally varied to explore a range of MW halo densities at the LMC’s recent 
perigalacticon, n p . A broad range of /3’s were also explored for the mid-core case, and a few miscellaneous parameter choices 
were also explored (black stars, some beyond the plotted range). Beneath these stars, the analytic models of Section [3] are also 
shown. 
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then slowly descends further into the MW halo, picking 
up speed and facing an increasingly stronger ram pres¬ 
sure, an effect compounded by the increasing CGM den¬ 
sity as the LMC travels towards pericenter of its orbit. 
The theoretical Rgg soon shrinks to become smaller than 
the LMC’s initial extent. At this point, from face-on pro¬ 
jections of the gaseous disk, we see material torn off the 
disk. The stripping is not entirely smooth, but rather 
episodic. The disk’s extent often lies a substantial frac¬ 
tion of a kiloparsec beyond Rgg, defying the expected 
stripping. Then, on a timescale of a few Myr, this disk 
material beyond Rgg fragments globally, and disinte¬ 
grates, beginning at the leading-edge. Any remaining 
large-scale asymmetries are evened out within one rota¬ 
tion period of the LMC disk. Material pulled beyond the 
disk plane is first “folded back” before flying off the disk. 
As discussed below, this often enhances the appearance, 
from our viewing perspective on Earth, of an asymmetric 
HI column distribution. From slices of gas density, we see 
that material unbound from the disk often lingers behind 
the LMC disk, shielded from the wind and still bound to 
the system, though at unobservably low column. 

Figure |6] provides mock observations from a subset 
of representative simulations compared to the observed 
HIPASS HI column map in the LOS point of view (LOS 
POV). For the simulations this LOS POV is a cartesian 
reduction (sum) performed with yt of physical density 
into surface density along equally spaced rays throughout 
the domain orientated along the LOS 5-axis, as defined 
in Table 3. The resulting 2D surface density map is cen¬ 
tered on the LMC kinematic center, with a north-vector 
aligned with the LOS y-axis. For the observed HI col¬ 
umn, we interpolated the data onto an orthographic pro¬ 
jection with astropy’s World Coordinate System (WCS) 
centered on the LMC kinematic center. In this projec¬ 
tion, x- and y -vectors are tangent to the lines of constant 
RA and DEC passing through this center, the former 
anti-aligned with RA so it increases left-to-right. Dis¬ 
tances are then mapped to kpc via the measured dis¬ 
tance to the LMC (see Table E} and by approximating 
the map as locally flat. Thus the LOS POV for both the 
simulations and observations roughly corresponds to 2D 
maps in the LOS frame defined in Section [O] The error 
introduced by this mapping at the observed Rt ~ 6 kpc 
is O(0 3 ) ~ (6 kpc/50 kpc) 3 « .2%. 

The central panel of Figure [6] shows a mock HI ob¬ 
servation of our fiducial simulation, M-1.2 at present 
day, « 50 Myr past the LMC’s recent perigalacti- 
con. The observed stellar d istribution is also overlaid 
(Ivan der Marcl fe Cionil 1200 If) as contours on the panel. 
To produce these plots, we tagged our simulation’s gas 
disk with a color field, to identify where disk mate¬ 
rial would end up in the simulation. We then elimi¬ 
nated any cells where this disk gas had fallen below .03 
cm~ 3 , i n density, as thi s material is likely highly ion¬ 
ized dRahmati et al.ll2013f) . We then used yt to produce 
a column density projection of the surviving HI along 
the LOS axis between the solar neighborhood and the 
LMC’s kinematic center. This vector is labeled 5 los in 
Table 3, and was found in the simulation frame, SIM, 
with the coordinate transformation matrices described 
in Section 14.21 A planar projection was used for sim¬ 
plicity, though the opening angle of the observed LMC 
disk is 9 ss tan _1 (5 kpc/50 kpc) = 5.7°, and so using 


non-converging column rays introduces errors ~ .1 kpc. 
We then cut out any region of the projection where the 
total HI column falls below 10 19 cm -2 . A similar cut 
was made to the observed HI distribution, shown in the 
bottom-right panel of Figure [G] The present direction of 
the headwind experienced by the LMC (in a frame where 
the LMC is at rest) is shown in the upper-left corner of 
the panel. 

Figure [6] also includes a black ellipse on each simula¬ 
tion panel to denote the present day truncation radius 
predicted by the GG model, Rgg- To provide a di¬ 
rect numerical comparison to this, we take the “leading 
quadrant” of the LMC disk, defined from 90° - 180° in 
the LOS projection frame, and search along rays equally 
spaced in angle for the last point where the HI column 
drops below 10 19 cm -3 . We then fit a ellipse to these 
points, by minimizing the sum of the square of the radial 
residuals to each point. This ellipse is shown in black 
with a white outline. Both the fitted and Gunn/Gott 
ellipses have orientations and axis ratios constrained to 
match a circle in the LMC frame as would be seen from 
the LOS POV. 

Analysis across all simulations indicated that this mea¬ 
sure of truncation radius is systematically smaller (and 
closer to the prediction of GG) than measurements taken 
from face-on projections of the LMC disk (as seen from 
a fictional observer directly above the LMC disk plane). 
This discrepancy vanishes if we use a higher column 
threshold, such as 10 2 ° cm -3 , providing evidence that 
the difference is caused by extraplanar, “folded back” 
material that is less visible from the LOS POV. Fur¬ 
ther, the systematic difference between this measurement 
taken from the LOS view and the face-on view is substan¬ 
tially wider than the scatter across the simulations. Thus 
this is an important feature of our simulations completely 
unexplained by the analytic GG model that will have a 
measurable impact on our inferred diffuse CGM density. 

Visual inspection of this gas beyond the predicted GG 
radius but within R t reveals a rippled structure to its 
column density, suggesting the gas has been destabilized 
by the wind, and is perhaps en route to becoming liber¬ 
ated from the disk. Visual inspection also reveals that 
the trailing (bottom-right) edge of the disk has gas well 
beyond the fitted truncation radius. This asymmetry in 
the HI morphology is both a product of viewing per¬ 
spective (folded back material we can see on the trailing 
edge and not the leading edge) and an actual discrep¬ 
ancy in the truncation radius on the leading and trailing 
edges, as seen in face-on projections of the disk. Bear in 
mind that the simulated stellar distribution’s iso-surface- 
density contours are themselves ellipses of the same axis 
ratio and orientation as those plotted in Figure [6] Thus 
we indeed see an offset in observable HI on the leading 
and trailing edges unexplained by the simple GG model. 
Notably, the discrepancy disappears for the highest den¬ 
sity runs shown here (e.g. Maxwind). 

An important point to make regarding these mock 
observations is the absence of high-column gas beyond 
the disk, particularly in the direction of the Magellanic 
Bridge (bottom-right). This stands in stark contrast to 
the observed HI profile, where a significant quantity of 
gas lies in the bridge region, which current models at¬ 
tribute to t he work of tidal f orces from the LMC-SMC 
interaction dBesla et al.ll201(t) . The absence of high col- 
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Figure 6. Present Day HI column distribution for 8 representative simulations compared against the observed HI column 
(bottom right panel). Also overlaid in grey contours is the observed stellar distribution from 2MASS Ivan der Marel fc Cionil 
(120011 1. The name of the run, which can be matched to Table m is displayed in the upper-right corner. The pericentric MW 
gaseous halo density, n p in each simulation is displayed in the bottom right of each panel in units of 10“ 4 cm -3 , with the 
corresponding peak ram pressure experienced at this point in the orbit displayed in the bottom right, also in cgs. These two 
quantities tend to increase from top to bottom in this hgure. The observed truncation radius, found from htting the leading edge 
of the disk (upper left quarter) where the HI column falls to 10 19 cm~ 2 , is displayed as black ellipse with a white outline, whose 
axis ratio is constrained to match the LMC disk orientation found in vdM02. The corresponding truncation radius predicted 
by Gunn & Gott is also plotted as a thinner black ellipse. The black arrow in the upper left of each panel indicates the 
direction of the headwind experienced by the LMC at present day. Finally, the white cross denotes the LMC kinematic center 
as traced by the HI 


umn gas beyond the disk plane is a feature present across 
our runs, along all sight-lines. 

The runs highlighted in Figure [6] were chosen to vi¬ 
sually illustrate some points that further analysis across 
time and across the broader suite of 42 simulations cor¬ 
roborates. First, we note the obvious, by considering 
the top-right, central, and bottom-left panels, i.e. runs 
M— 0.8, M—1.2 and M— 3.0, which feature the same cen¬ 
tral MW gaseous halo density but decreasing /3, and by 
extension larger n p and higher ram pressure: consistent 
with the model of GG, these runs feature progressively 
smaller truncation radii. 

Next, we note the central row shows runs H-l. 2, M-l. 2 
(fiducial), and L-l. 2, which all exhibit identical pericen¬ 
tric densities but decreasing /3 (i.e., from left to right 
these MW gas halos have slower density falloffs). From 
visual inspection, the runs appear more or less identical, 
and the discrepancy in their fitted truncation radii do 
not follow any trend with /3. These anecdotes suggest a 
point we strengthen later: that the LMC’s present day 
HI truncation radius is a clean probe of local density, 
n p = n(R ~ 48 kpc), i.e. the maximal ram pressure sets 
the value, and is independent of the broader halo profile 
the system has orbited through. These simulations thus 
bolster a key assumption of Section [3] 

Another comparison worth highlighting involves the 


top-left panel and the central fiducial panel. The for¬ 
mer shows a face-on wind simulation, FM-1.2, where we 
launched the wind directly along the z-axis of the LMC 
disk plane, eliminating any asymmetry between how dif¬ 
ferent quadrants of the disk plane experience the head¬ 
wind. This run is otherwise identical (in wind speed, 
density, temperature and evolution) to the fiducial M-l. 2 
run. The first important point to make is the inferred 
truncation radius of 6.91 kpc is substantially larger than 
any of the values in the central row of Figure[6l despite all 
these runs having identical n p , putting it in even further 
disagreement with the GG model’s prediction. The sec¬ 
ond important point involves asymmetry in the observed 
HI column: despite the symmetric pummeling of all edges 
of the disk, from the LOS POV the gas still appears sub¬ 
stantially farther from the LMC kinematic center on the 
trailing edge of the disk. This is the clearest example yet 
of extra-planar, “folded back” material influencing the 
observed HI distribution. It further highlights the im¬ 
portance of viewing perspective. From visual inspection 
we see this gas is at substantially lower column (though 
still above our cut) and exhibits a wavy, unstable struc¬ 
ture rife with filaments slowly breaking away from the 
disk. 

Finally, we call attention to the top-center panel, run 
M-0.13 where the typical relation between Rgg and Rt 
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has been reversed: the simulated truncation of the mock 
HI observation lies within the GG model’s prediction. 
This is simply a product of our density cut choice: be¬ 
fore the wind has even gained enough strength to induce 
any stripping of LMC disk material, gas below our phys¬ 
ical density cut of .03 cm -3 dominates at large radii. 
Whatever high density gas does go into our mock den¬ 
sity projection is then removed by the surface density 
cut at 10 19 cm -3 . Indeed, choosing a cut of 10 18 restores 
agreement, though the probed radius exhibits an even 
higher degree of scatter across the simulations and dis¬ 
agrees even further with the prediction of the GG model. 
This issue motivates our choice of a higher column cutoff. 
However, if we choose too high a surface density cut — 
say, 10 2 ° cm -2 — the observed truncation radius of the 
initial HI mock observation is smaller than the present 
observed extent of the LMC, rendering our predictions 
useless. Thus 10 19 cm -2 is optimal for our purposes. 
Varying the physical density cut from .03 cm -3 likewise 
impacts the results, though this value is motivated by 
studies of the ISM and not an arbitrary choice. 

4.5.2. Radial Profiles of HI Column 

Figure [7] quantifies the surface density profiles of Fig¬ 
ure [G] by sampling the mock-HI surface density along a 
dozen rays from the LMC disk center out towards the 
edge in each disk quadrant (clockwise from top-left: the 
leading edge, southeast, towards the bridge/stream and 
northwest). The 2D distance along these rays is then 
mapped to a radial distance from the LMC kinematic 
center assuming all material lies in the disk plane, and 
then these profiles are added together into an average 
surface density profile along the leading edge. We per¬ 
formed this analysis in an identical fashion for both the 
simulations and the observed LMC HI distribution. From 
this we find the runs with n p = 1.2 x 10” 4 cm” 3 (regard¬ 
less of hi/mid/lo core model) produce the most successful 
match to the observed surface density profile along the 
leading edge. The agreement among these three simula¬ 
tions further bolsters our argument that RPS is a clean 
probe of n p , rather than the broader halo profile. The 
face-on wind model at the same pericentric density mean¬ 
while produces consistently higher HI surface density, 
consistent with Figure [6] Also shown is a substantially 
higher density MaxWind model, which truncates the disk 
at much smaller radius. 

We should emphasize here that our goal was to match 
the observed truncation radius, and not the broader HI 
profile. For the former goal, we considered where the col¬ 
umn dropped to 10” 19 cm” 2 , a probe of compression and 
oblation of gas by ram pressure stripping. The broader 
profile however ought to be sensitive to the details of 
compression-induced star formation, as HI may be pro¬ 
cessed into stars, outflows and ionized gas, dropping the 
observed column deeper into the disk. This could per¬ 
haps even account for the steeper falloff of the leading 
edge, beginning at r ~ 4 kpc, which our simulations do 
a poor job recreating. 

Beyond the leading edge, we find the simulated gas 
disk extends further, hitting the column threshold of 
10” 19 cm” 2 beyond 7 kpc in the remaining three quad¬ 
rants. Comparing to the observed profile, the southeast 
quadrant (bottom-left) features the closest match, ignor¬ 
ing the 30-Doradus anomaly in the inner profile. In the 


southeast quadrant leading in to the MS and bridge, the 
simulated profile fails to reproduce the high column at 
large radii seen in the observed system, which is expected 
since interactions with the SMC are not simulated. 

4.5.3. Time Evolution of Leading Edge Truncation Radius 

We next explored the time-evolving nature of the trun¬ 
cation radius, found from mock LOS HI column maps 
such as in Figure [6l by again fitting the radius at which 
the gas column density drops to 10 19 cm” 2 and then 
fitting an ellipse whose shape comes from assuming cir¬ 
cular contours in the LMC disk plane. Figure [5] shows 
how this fitted truncation radius evolves over the final gi- 
gayear of the LMC’s orbit, across a suite of simulations, 
and compares this to the radius predicted by the GG 
model assuming instantaneous stripping. 

The top-left panel of Figure [8] plots the fiducial simu¬ 
lation along with a host of higher and lower density halo 
models, all from the “mid-core” set (i.e. MW gas halos 
with the same core radius and density, but different ex¬ 
ponential falloffs). All runs exhibit oscillations at early 
times, with amplitudes ~ 10% of the HI truncation disk 
radius. These oscillations are leftover from our initial 
conditions, and are similar to those found in runs devoid 
of radiative cooling in T09. When the wind becomes 
sufficiently strong as to alter the observed truncation ra¬ 
dius, there is an initial rise in the truncation radius. This 
is due to vertical compression of the disk, which causes 
more gas to lie above our p = .03 cm” 3 density cut 
since similar plots that did not make such a cut do not 
exhibit the same spike. From there, the truncation ra¬ 
dius shrinks rapidly as the LMC approaches pericenter, 
its slope in the radius-time graph consistent with the an¬ 
alytic GG prediction, though its observed value is about 
5% larger than the model predicts for mid-density runs 
(such as the fiducial case). Some of the runs, the fidu¬ 
cial run among them, then pull up further from the GG 
prediction near present day, which we discuss more while 
describing the next panel. 

The bottom-left panel of Figure [8] shows six runs with 
an identical MW halo density at pericentric passage. An 
instantaneous application of the GG model at pericenter 
would thus predict an identical truncation radius across 
these models. Indeed runs with three different central 
halo densities and disparate /3s — H-l. 2, Fiducial, and 
L-1.2, in red, black and blue, respectively — do not 
show a significant difference in the inferred present-day 
truncation radius. This result illustrates that the LMC’s 
present day truncation radius is indeed a localized probe 
of the MW gaseous halo’s density. 

Figure [8js bottom-left panel also plots hi-, mid- and 
lo- core runs at the fiducial n p for simulations with a 
face-on wind. While the early evolution for all three core 
models agrees well, the inferred truncation radius departs 
sharply from the full-orbit winds at late times, with the 
truncation radius rising rapidly in two of the three face- 
on wind runs. Although we do not explore the cause of 
this disparity in detail, it seems the symmetric nature of 
the wind perpendicular to the galaxy plane has sustained 
and enhanced, rather than damped-down oscillations in 
the inferred HI truncation radius. 

The final, right panel of Figure [H] plots the ratio of the 
simulated truncation radius to that predicted by GG, 
to quantify the disparity. Blue and gold shaded regions 
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RADIUS [ RPC ] 


Figure 7. Averaged profile of surface density in the four quadrants of the LMC, computed as in Figure Q] Black lines shows 
our fiducial simulation, whereas red and blue lines show hi- and lo-core halo density models with an identical pericentric density 
of 1.2 x 10” 4 cm -3 . A thick grey line denotes the same average profile found from the observed HI column map. Also shown 
are the higher-density MaxWind simulation (purple) and a face-on wind run, FM-1.2 (orange). 


show the bounds within which all of our simulations lie, 
across various core concentration and pericentric density 
choices, for the full-orbit and face-on winds respectively. 
From this we find that the full-orbit simulations produce 
a truncation radius ~ 5% higher than GG across all runs, 
with a comparatively small scatter, whereas the face-on 
wind leads to radii over 10% higher than the analytic 
model. These oscillations have a period of roughly 200 
Myr, whereas the period of epicyclic oscillations in the 
disk is « 450 Myr at 6 kpc. 

4.5.4. Comparisons between GG toy model and simulations 

To move towards a quantitative prediction for the 
MW’s diffuse CGM density, we next sought to quantify 
the inferred HI truncation radius across our entire suite 
of simulations. For each run, we produced the mock HI 
observations, as in Figure [G] and once again fit a con¬ 
strained ellipse to the column map where the surface 
density dropped to 10 19 cm” 2 . Figure [9] shows the re¬ 
sults of this effort, plotting inferred truncation radius 
against each orbit model’s pericentric MW halo density, 
n p . From this we found a relation whose trend was well 
matched to the analytic prediction of GG, though near 
the observed truncation radius the distribution was con¬ 
sistently offset by ~ 5%, just larger than the scatter in 


the results. Our “mid-core” halo models sampled roughly 
two orders of magnitude in n p , with a dense sampling 
around the 1.2 x 10” 4 cm” 3 fiducial run, where the re¬ 
sults best matched the observed truncation radius. The 
simulations break from the GG model at low density, 
since for our mock HI observations the density cut at 
.03 cm” 3 truncates even our initial disk near 7 kpc. In 
addition to sampling across n p , we ran numerous sim¬ 
ulations with more concentrated core densities (hi-core, 
red stars) and less concentrated (lo-core, blue stars). Al¬ 
though small trends are evident within and among these 
three model classes, they exist within the general scat¬ 
ter, suggesting the LMC HI truncation radius along the 
leading edge is indeed insensitive to the broader halo gas 
profile, and is thus a clean probe of the gas density at 
ss 48 kpc. This conclusion sta nds i n stark contras t to the 
orbital scenario explored by iMastropietro et all (I2009T ) 
whose lower-speed orbital scenario suggested the strip¬ 
ping process probed the broader gas profile over a longer 
duration of the LMC’s orbit. A pair of our exploratory 
runs (T-6 and T-6X, purple stars) really serve to em¬ 
phasize this locality result: these two models were both 
chosen to produce an n p that would truncate the disk 
at exactly 6 kpc, as predicted by the GG toy model. 
But they were chosen from opposite ends of the /3-model 
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Figure 8. The top left panel shows the evolution over time of the leading edge truncation radius, found from LOS column maps, 
in our simulations. The thick black lines denotes our fiducial simulation, with a pericentric density n p = 1.2 x 10~ 4 cm~ 3 , while 
green solid lines show a variety of other pericentric densities. Dashed lines show the instantaneous stripping radius predicted 
by the GG model. The bottom-left panel shows the same format as the top plot, but for models with identical pericentric halo 
densities but different exponential falloffs. The red curves are for the hi core, concentrated model, whereas the blue curves are for 
lo core, extended model. Also shown as dash-dot lines are so-called “face-on” runs, where the wind was launched perpendicular 
to the galaxy’s disk plane, rather than as it would be experienced by an orbiting LMC. Finally the right plot shows the ratio 
of the simulated truncation radius to that predicted by GG. Blue and gold swatches denote the range spanned by all full-orbit- 
and face-on-wind runs, with the fiducial halo models highlighted in black. 


parameter space. Despite this, they appear nearly over¬ 
lapping in Figure[9j with a separation much smaller than 
the overall scatter across our runs. 

The right panel of Figure [8] suggests that the wind’s 
angle of approach (which we denote ip for brevity) might 
have a measurable impact on the implied MW halo den¬ 
sity. Flere ip(t) is the angle between disk plane’s normal 
vector and the wind velocity vector. In the majority of 
our simulations this quantity varies over time as the LMC 
orbits the MW. We assume the angular momentum vec¬ 
tor is stable for the fina l 200 Myr of the orbit (consistent 
with iBesla ct all 12012ft , when the ram pressure ramps 
up to where it measurably impacts the disk gas. Under 
this scenario ip begins at 20°, climbs past 40° at peri- 
center, and finishes off at the observed present value of 
about 50°. For our face-on simulations, by construction, 
ip it) = 0, even as the wind’s speed and density evolve 
as in the full-orbit runs. These face-on runs showed the 
ip = 0 runs have a systematically higher offset from the 
GG model, due to larger oscillations in the disk radius 
post-stripping than the full-orbit ip p ~ 40° runs. This 


suggests the disk’s instantaneous truncation radius at 
present day may be sensitive to this angle. Thus we 
were led to perform a few simulations with a constant 
intermediate wind orientation of ip = 30°. We ran three 
simulations, all at the fiducial n p = 1.2 cm -3 , across our 
three core models. These inclined wind runs are desig¬ 
nated IL-1.2, IM-1.2 and IH-1.2. Figure flOl plots these 
runs against both the full-orbit and face-on simulations. 
We see that the 30° inclined wind runs across the three 
core models agree almost precisely with their full-orbit 
counterparts. This suggests the symmetry of the face-on 
run leads to enhanced oscillatory behavior qualitatively 
distinct from runs at any intermediate inclinations, and 
that our results remain insensitive to the exact value of 
the wind angle, ip. 

4.5.5. Influence of Ram Pressure Stripping on Velocity Field 

of HI disk 

Beyond understanding how RPS alters the spatial dis¬ 
tribution of the LMC’s HI, our simulations allow us to 
study the gas’s velocity structure. Figure [lT] again pro¬ 
vides mock HI observations, this time of the radial (LOS 
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Figure 9. The leading edge truncation radius, as inferred from mock HI observations such as in Figure [6j as a function of 
each simulation’s pericentric MW gaseous halo density, n p , for our entire suite of simulations using our standard LMC model 
described in Section o Stars denote each simulation, with orbit models spanning two orders of magnitude in n p , with a denser 
sampling where the simulation’s results well matched the observed leading edge truncation radius (horizontal dashed line, with 
grey region denoting a 3 — er confidence region). In addition to sampling across pericentric densities, we also explored whether 
or not the broader halo profile impacted the observed radius, with the lo-, mid-, and hi-core models, displayed here in blue, 
green and red respectively (see Section [4.411 . Although trends within and among these three classes of halo model are evident 
here, they exist within the general scatter, suggesting the observed truncation radius is indeed a clean probe of the local MW 
gaseous halo density at r = 48 kpc from the Galactic center. 


component of) ve locity, as well a s the observed HI ra¬ 
dial velocity from [Putman ct Till ( 2003 ). To form these 
images, we again sought to eliminate gas that is likely 
ionized, removing gas where n < .03 cm -3 . For each 
pixel we took a density-weighted average of uloSj the 
component of the gas’s velocity towards or away from 
the observer. Finally, we again eliminated from the map 

any region where the gas surface density fell below 10 19 
_2 
cm . 

For the observed HI, the raw ulos data needed to be 
corrected for both the LMC’s proper motion through 
space (relative to the Sun) and the disk plane’s solid 
body motion about its center of mass. Only then could 
a fair comparison be made to our simulated disk which 
sits at rest in the “wind tunnel”. To achieve this, we 
subtracted vel ocity corrections using m achinery laid out 
in vdM02 and Ivan der Marel fe Cionil (I2001D , described 
in Appendix lAl 

For reference, Figure [Tl] also plots the observed stel¬ 
lar distribution, the truncation radius found as in Fig¬ 
ure E and the GG prediction for the truncation radius. 
Three simulations are shown: our fiducial (M-l. 2) run, a 
face-on wind run with an identical halo density model to 
the fiducial (FM-1.2) and our highest density wind run, 


MaxWind. For all three runs, regardless of wind direction 
or strength, the ram pressure fails to strongly perturb 
the basic gradient of a circularized, rotating disk. In con¬ 
trast, the observed radial velocity in the disk plane is rife 
with substructure, most notably a swirl pattern with sev¬ 
eral arm-like structures. Because this observed velocity 
map’s appearance depends on accurate parameters for 
the LMC’s three-dimensional position, orientation and 
bulk motion, we varied these parameters systematically 
and found no choice of parameters that could severely 
diminish these swirls. We conclude the LMC disk plane 
indeed features a rich velocity structure unexplained by 
its proper motion, solid body rotation or circular mo¬ 
tion. Our simulations also suggest that RPS cannot be 
responsible for this irregular motion to any appreciable 
extent. 

4.6. Implications of an LMC disk with a Higher Initial 
Gas Fraction 

In this section we consider alternative models for the 
LMC’s gaseous disk profile, and explore how they inform 
the implied CGM density. The LMC gas disk model 
chosen in Section [2. II represents perhaps a lower-limit on 
the LMC disk’s gas fraction prior to its recent pericen- 
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Gas Rich Disk Simulations 


n 0 

r c 

P 

p P 

Total Mass 

n p 

S GG 

R'j 1 

Name 

0.139 

0.100 

0.323 

4.07e-13 

6.46e+10 

3.50e-04 

6.63 

7.32 

GL-3.5 

0.139 

0.100 

0.316 

4.65e-13 

7.66+10 

4.06e-04 

6.42 

6.51 

GL-4.0 

0.139 

0.100 

0.309 

5.23e-13 

8.71 + 10 

4.54e-04 

6.21 

6.76 

GL-4.5 

0.139 

0.100 

0.278 

9.29e-13 

1.69e+ll 

8.06e-04 

5.32 

5.75 

GL-8.0 

0.139 

0.100 

0.275 

9.88e-13 

1.80e+ll 

8.52e-04 

5.23 

5.36 

GL-8.5 

0.139 

0.100 

0.272 

1.05e-12 

1.92e+ll 

9.01e-04 

5.15 

5.21 

GL-9.0 

0.46 

0.35 

0.559 

1.39e-13 

1.25e+10 

1.21e-04 

8.39 

8.73 

GM-1.2 

0.46 

0.35 

0.486 

4.07e-13 

4.25e+10 

3.51e-04 

6.62 

7.35 

GM-3.5 

0.46 

0.35 

0.477 

4.65e-13 

4.96e+10 

4.01e-04 

6.41 

7.01 

GM-4.0 

0.46 

0.35 

0.469 

5.23e-13 

5.69e+10 

4.52e-04 

6.22 

6.83 

GM-4.5 

0.46 

0.35 

0.430 

9.29e-13 

l.lle+11 

8.03e-04 

5.32 

5.88 

GM-8.0 

0.46 

0.35 

0.426 

9.88e-13 

1.19e+ll 

8.52e-04 

5.23 

5.99 

GM-8.5 

0.46 

0.35 

0.422 

1.05e-12 

1.28e+ll 

9.04e-04 

5.14 

5.88 

GM-9.0 

0.810 

0.800 

0.630 

4.07e-13 

3.18e+10 

3.52e-04 

6.62 

7.10 

GH-3.5 

0.810 

0.800 

0.620 

4.65e-13 

3.71e+10 

4.03e-04 

6.42 

6.96 

GH-4.0 

0.810 

0.800 

0.610 

5.23e-13 

4.27e+10 

4.55e-04 

6.23 

6.91 

GH-4.5 

0.810 

0.800 

0.563 

9.29e-13 

8.19e+10 

8.02e-04 

5.32 

5.96 

GH-8.0 

0.810 

0.800 

0.558 

9.88e-13 

8.80e+10 

8.52e-04 

5.23 

5.87 

GH-8.5 

0.810 

0.800 

0.554 

1.05e-12 

9.46e+10 

9.06e-04 

5.16 

5.80 

GH-9.0 


Table 5: Our 19 gas-rich LMC disk model simulations (all prefixed with G). These runs are grouped by MW gaseous halo core 
density as in Table l4~4l ( GL-, GM- and GH- runs) with various /3 chosen to explore a range of pericentric densities. See Table 
521 for a description of column quantities. 
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Figure 10. A comparison of runs by the angle between the 
wind and disk plane’s normal vector, ip, plotting quantities as 
in Figure [9] Black stars denote runs where the wind velocity 
mimics the orbit of the LMC, with ip(t) ~ 20° —¥ 50°, with 
a pericentric value ~ 40°. Yellow circles denote face-on runs, 
with ip = 0 at all times, even as the wind speed and density 
vary as in the full-orbit runs. Finally, triangles show three 
runs at an intermediate, fixed ip = 30°. Blue, green and red 
markers denote IL-1.2, IM-1.2, and IH-1.2, for our three core 
density models with n v = 1.2 x 10~ 4 cm -3 . While we find 
the face-on wind produces consistently higher inferred trun¬ 
cation radii at present day (due to the oscillations observed 
in Figure EJ, both the 30-degree and full-orbit (ip p « 40°) 
runs show excellent agreement in R t , well within the scatter 
of our simulations. Thus we conclude the exact nature of disk 
orientation is of minimal importance, so long as the wind is 
not face-on. 

trie passage. Alternative models that still obey the mass 
constraint of K98 are possible, producing gaseous and 
stellar disk properties (e.g. extent relations and gas frac¬ 
tions) consisten t with observations of isolated late-type 
dwarf systems (ISwaters et al.l 120021 : iBeeum et all 120051 : 
iKreckel et al.ll201lh . As the gas mass increases, the pro¬ 
files do a consistently worse job matching the LMC’s 
present-day observed surface density profile, producing 
columns too low at its center with too shallow a radial 
falloff. However, ram pressure dynamics, including con¬ 
tour compression and enhanced SF can perhaps alleviate 
these issues at late times, especially since the headwind 
required to truncate the HI disk at its observed extent 
grows increasingly strong for these alternative models. 
The model we consider boosts the initial LMC gas mass 
by over a factor of two, to ~ 10 9 M 0 . Recent UV absorp¬ 
tion line data coupled with Cloudy modeling of ionized 
gas in the MS suggest the amount of gas strip ped from 
the M Cs may exceed their current HI mass (|Fox et al.l 
mu, so this scenario is not outlandish. 

To explore the scenario of an Hi-rich initial LMC, we 
formulated a second set of initial conditions with an 
identical stellar and DM distribution, but with a gas 
disk whose scale radius was twice as extended as the 


standard run’s scale radius (and thus twice the stars): 
Ogas = 3.4 kpc. We then set the mass normalization con¬ 
stant, M gas = 1.22 x 10 9 Mq, which once-again ensured 
the dynamical mass within R = 4 kpc matched the con- 
stra i nt of K98 . This result also fell on the relation of 
ISwaters et al.l (120021b suggesting it as a plausible scenario 
for a gas-rich late-type dwarf galaxy before plunging into 
the MW halo and undergoing RPS. This produced an 
initial mock-HI surface density that was too low at the 
galaxy’s center and much too high towards its edges in 
all directions, but this in itself does not preclude such a 
model, since the wind required to strip such an LMC to 
its present observed extent will be much higher density, 
and thus has the potential to dramatically reshape the 
entire HI profile. 

We again performed a suite of simulations with MW 
gaseous halo models spanning a wide range of n p , and 
also varied the halo’s core density, as before. Table 
5 summarizes these runs. Our n p choices are concen¬ 
trated around where the predicted truncation radius 
would roughly match the present-day observed value ~ 6 
kpc, but a run was also done at the fiducial density of 
n p = 1.2 x 10 -4 cm -3 for comparison. These runs are 
prefixed GL-, GM- and GH-, with the first letter denoting 
their “gas rich” initial disk and the second their halo core 
model. 

Figure fl2l shows our canonical mock-HI observation as 
before for a representative gas-rich run, GM-8.5, which 
produced a stripping radius in rough agreement with the 
standard disk’s fiducial run and the observed distribu¬ 
tion (also repeated in the new figure for comparison). 
Immediately evident is the higher gas column towards 
the disk edge, compared to the standard simulation and 
the observed map. This is seen most easily in the bot¬ 
tom panel of Figure [T2l where we again plot the averaged 
radial profile along the leading edge (as in Figure 0. In 
this sense, we find RPS cannot substantially alter the 
global gas profile, even when the wind has been ramped 
up to sufficient density to match the leading edge’s ex¬ 
tent. The rate of star formation required to alleviate this 
discrepancy is ~ lOMg/year, unrealistically high for the 
LMC, given observations. From visual inspection of the 
column maps that the gas-rich run’s flocculent disk edge, 
shaped by the denser, more violent halo headwind, does 
a better job matching the observed profile. But here we 
caution that a more detailed treatment of ISM physics 
could easily reverse these fortunes. In all, the gas-rich 
runs provide a poor match to the observed HI content of 
today’s LMC, and thus we deem this such high gas disk 
models to be unlikely. 

We can once again move towards an inferred halo den¬ 
sity and compare to the standard runs by considering 
the whole suite of gas-rich simulations. Figure [13] ac¬ 
complishes this, plotting the same derived quantities as 
in Figure [9] This figure also compares both gas-rich and 
standard disk runs to the GG model’s analytic predic¬ 
tion. The higher surface density gas disk shifts the GG 
model substantially, to give an inferred MW gaseous halo 
density ~ 5 times higher than the runs with the standard 
disk ICs. The simulated results follow this trend, once 
again exhibiting an offset from the analytic model and 
a scatter roughly equivalent to the original simulation 
suite. Figure [13] also plots the ratio of the simulated 
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Figure 11 . Comparison of simulated HI radial velocity from mock-HI observations to the observed HI radial velocity (see text 
for explanation of proper motion corrections). Red and blue indicate motion away from and towards the observer, whereas 
white indicates gas with no radial velocity. Yellow contours again over plot the observed stellar profile llvan der Marel fe Cionil 
1200111 . and white and dashed-line ellipses denote the truncation radius fitted to the HI column profile and predicted by the GG 
model, respectively. The simulations shown here are for the face-on and fiducial runs with n p = 1.2 x 10 -4 cm' 3 (top-left and 
bottom-left) and the “max wind” scenario with n v = 8.44 x 10 -4 cm' 3 (top-right). 


truncation radius, Rt , to the GG model, Rgg- Here we 
see the gas-rich runs feature an offset from GG of ss 5% 
on average, in good agreement with the standard runs. 
Note that this GG model includes the full self-gravity of 
the disk gas, typically ignored in RPS models of MW- 
sized systems with lower gas fractions. Including this 
self-gravity was captured by setting a = 1.0 in Equation 
03 Other choices for a, in particular the standard choice 
of a = 0, led to larger offsets between the GG model and 
the simulations, most for the gas-rich runs. This sug¬ 
gests including the full self-gravity of the gas is the most 
accurate approach. 

5. INFERRING THE MW GASEOUS HALO DENSITY AT 
LMC DISTANCES 

The purpose of this section is to extend the analytic 
considerations of Section [3] with the simulation results of 
Section 0] to produce a MW CGM density bounds im¬ 
plied from the truncation radius, Rt, of the LMC HI 
disk’s leading edge. To do so, we first characterize the 
offset between the CGM density implied by the analytic 
GG model, and that found from simulations, as cap¬ 
tured by the offset parameter D in Equation [10] This 
together with a systematic enumeration of model uncer¬ 
tainties permits us to produce a quantitative prediction 
for n p (48 kpc), the density of the MW diffuse gaseous 


halo at the distance of the LMC’s most recent pericen- 
ter. We begin by describing the results of this effort, 
before going back to enumerate the various uncertainties 
we chose to model. 

5.1. A Density Constraint for the MW CGM at r = 48 

kpc 

Using Equation [TU] together with samples from the er¬ 
ror space of its various parameters (described below) we 
obtained an implied MW CGM density of 

n M WHaio (R = 48.2 ± 2.5 kpc) = l.lt ^ x 1(T 4 cm" 3 . 

(13) 

Figure [14] displays the results of our full model graph¬ 
ically. The upper panels display the spread in the sam¬ 
pled Rt and computed values n p , P p and MMWgashaio- 
The first quantity was directly sampled from a normal 
distribution R t ~ N( 6.2 kpc, (.25 kpc) 2 ). To judge the 
importance of this choice, we also used a uniform distri¬ 
bution with an identical standard deviation, which had 
only a small effect on the derived shapes and spreads in 
our other quantities. The pericentric MW gaseous halo 
density, n p (and corresponding ram pressures) found by 
our full model are substantially higher than the results 
of Section [3] a reflection of the offset we found in our 
simulations, thus implying a more gas-rich MW halo at 












24 


Salem et al. 




RADIUS [ KPC ] 


Figure 12. Results for our gas-rich LMC disk model simulations, which began with a higher initial gas mass of M gas , tot « 
1.8 x 10 9 Mg. Top: mock HI column map as in Figure [9]for the new gas-rich fiducial halo model run, GM-1.2 (center panel), 
compared to the fiducial simulation, M-1.2 (left), and observed HI distribution (right). Producing f?t ~ 6 kpc now requires a 
MW diffuse halo density at 50 kpc a full factor of 8 higher than for our standard disk model. This more violent scenario leads 
to a splotchier disk edge, which appears to match the observed profile better. However, ram pressure fails to alter the the gas 
column at intermediate radii, which remains noticeably higher than the observed system’s profile, as the bottom panel’s plot 
of the leading edge surface density profile (as in Figure Q shows. This panel also plots the fiducial simulation (black line) which 
features a falloff more in line with the observed system. These gas-rich results would require an unrealistically high amount 
of SF and feedback to bring this profile in line with current observations, suggesting such a pre-infall LMC HI disk is rather 
far-fetched. 


~ 50 kpc than the toy model suggests, by a factor of 
two. We can go a step further, albeit with healthy dose 
of extrapolation. Section 16.11 extrapolates this density 
result to compute a total mass for the MW CGM, and 
compares to recent models based on studies in emission 
and absorption. 

The final lower panel of Figure Q31 shows the agree¬ 
ment between our full model results and those of MB 15 
in parameter space of the /3-model for the MW halo. 
From this we find, across halo core density, our model 
predicts a consistently shallower falloff in halo density 
(smaller /?), implying a denser CGM profile. This de¬ 
gree of disagreement is hardly surprising: the hydrody¬ 
namic interaction of RPS involves the entire multiphase 
CGM, not merely the hot, X-ray emitting component. 
Fin ally, our R PS studies, as well as previous work (e.g. 
iRoedieer fc Henslerl 120051) shows that (supersonic) RPS 
is insensitive to the temperature of the ambient CGM, 
whereas X-ray models rely on thermodynamic assump¬ 
tions for the gas. However, given these considerations, 
the degree of overlap we do find with MB15 suggests 
both probes are converging on an accurate picture for 
the MW’s diffuse halo gas at intermediate radii. 


5.2. Sampling the Offset and Model Uncertainties 

To produce the above density bound, we needed to cal¬ 
ibrate the offset term D in Equation [TO] and then sample 
D along with a host of other uncertain parameters. We 
finish this section by enumerating the sample distribu¬ 
tions for these parameters, appealing to our full suite of 
simulations presented in Section 14.51 and observed prop¬ 
erties of the LMC. 

Simulation Scatter (D ):— The scatter observed in the 
simulated disk truncation radius is primarily due to os¬ 
cillations in the disk extent, as seen in Figure [51 which 
we choose to model as simply random noise, as Figure 
[5] indeed shows the phase of the oscillations shifts con¬ 
siderably across runs with no clear trend. Our data 
suggested the offsets in logspace did not have a strong 
trend with truncation radius (a t-test for a correlation 
coefficient produced a p-value of 33%) and are consis¬ 
tent with being drawn from a normal distribution (a 
Shapiro-Wilkes test yielded a p-value of 13%). We thus 
assume {D} ~ N(B,aff) 1 where B is the mean of the 
offset and <j 2 d its variance. We implicitly assume the 
offsets are independent of other changes in the LMC or¬ 
bit and disk structure parameters. We then used emcee 
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Figure 13. A comparison of our gas-rich initial LMC disk simulations (squares) to our standard disk model simulations (stars). 
The left panel plots HI truncation radius, Rt, versus pericentric halo gas density, n v , as in Figure[9] The GG model is shown 
now for both sets of ICs (two black lines). The observed HI truncation radius is also shown as a dashed line. The systematic 
offset seen in the standard runs is again observed in the gas-rich set, with an equivalent offset in n p . From the simulations we 
find the gas-rich model requires a MW halo gas density « 7 times higher than the standard runs at ~ 50 kpc. The right panel 
takes a ratio of Rt to the GG prediction across all simulations, with unity indicating a perfect match (dotted line). The gas rich 
runs feature an offset statistically consistent with the original runs. This agreement is strong evidence for the choice of a = 1.0 
in Equation [9l since models with a diminished gas self-gravity (a < 1.0) provided a worse match with disparate offsets between 
gas-rich and standard runs. 


(jForeman-Mackev et al.ll2013f) to find the joint distribu¬ 
tion of B and ctd, from which to sample moving forward. 
This model is a quick way to capture uncertainties in the 
simulations that include everything from the aforemen¬ 
tioned oscillations to the crude nature of mapping a 3D 
turbulent gas structure to a single truncation radius at 
a snapshot in time. 

Orbital Uncertainties :— Our work has relied on a fidu¬ 
cial orbit taken from a large sample of backwards-orbit 
integration calculations for the LMC/MW system that 
makes use of the latest proper motion measurements 
of K13. The integrations treat the LMC as a Plum¬ 
mer profile within an NFW profile for the MW halo, 
completely ignoring the SMC’s gravitational influence; 
these assumptions do not appreciably degrade the calcu¬ 
lation’s integrity during the past ^100 Myr during which 
the LMC undergoes its most recent peric entric passage 
(IKaHivayaiil et all 12011: iResla, et alll2007lh 10,000 MC 
drawings were made for each combination of Mmw G 
{1,1.5, 2 } x 10 12 M 0 and M LM c G {.3, .5, . 8 ,1,5,8} x 
lO n M 0 , though the inclusion of multiple galaxy masses 
enhances the spread in orbital parameters by less than 
5%. From these samples we find r p = 48.1 ± 2.5 kpc, 
v p = 340 ± 19 km/s and t p = 46.4 ± 8.5 Myr, where 
t p indicates time since pericentric passage. Our numer¬ 
ous simulations show the observed truncation radius is 
roughly independent of the broader halo gas density pro¬ 
file, and thus a well localized probe of pericentric halo 
conditions. Thus the uncertainty in r p , the distance of 
perigalacticon from Sgr A*, will not impact our den¬ 
sity measurement but rather directly translates to uncer¬ 
tainty in where within the MW halo that measurement 


was made. In contrast, uncertainty in v p is very impor¬ 
tant to n p , since the GG model depends quadratically 
on velocity. We thus sample this value from the orbit 
calculations when constructing our final measurement, 
along with the simulation scatter. To judge the impor¬ 
tance of the timing errors, we took measurements of the 
simulated truncation radius both 20 Myr beyond and be¬ 
fore the present day of our simulations. This introduced 
negligible scatter, and thus we ignore timing errors. 

Observed Truncation Radius (Rt ):— Section [3] describes 
how we inferred our measured Rt = 6.2 ± .25 kpc. To 
reiterate, our measurement of the radius of the LMC’s 
observed leading edge ignores the flocculent nature of the 
HI profile, punctured by holes caused by star formation, 
accretion flows and other phenomena of the turbulent 
ISM. We thus sample Rt from a normal distribution with 
a = .25 kpc, a conservative (large) estimate for accuracy 
of the leading edge’s position. As previously stated, we 
have neglected radiative cooling and star formation in 
this analysis as T09 found agreement within 10% of GG 
when accounting for a multiphase ISM. We further note 
that supergiant holes exist in the ISM of the LMC with 
sizes of on average 0.5 kpc (iKim et al.ll2003l) . This size 
scale is encompassed within the uncertainty of the trun¬ 
cation radius of the LMC outer disk. In this analysis 
we have also assumed that the MW CGM is smooth. If 
instead the LMC were primarily stripped by a clumpier 
medium, the clumps would need to be ~20 kpc in size in 
order to explain the roughly uniform truncation of the 
LMC’s leading edge (which had an initial disk radius of 
~10 kpc to match the current distribution of stars). En¬ 
counters with smaller, dense clumps may lead to a similar 
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Simulation / Full Model Results 
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Figure 14. TOP: Distribution of important LMC RPS model params from our full model incorporating uncertainties described 
in Section [5] including offset from GG model found from our suite of simulations. Left to right and top to bottom, these are 
the observed HI truncation radius (modeled as a gaussian, though this choice does not strongly influence results); MW gaseous 
halo density at LMC’s pericenter, n p = n(r p = 48.2 ± 2.5 kpc); and ram pressure experienced by the LMC at pericenter. To 
compute this final value, we also sampled uniformly from ngr^ G [.005, .03], motivated by MB15’s error space. The average 
of this inferred halo mass is a factor of three higher than the value found in the toy model of Section [3] BOTTOM: Colored 
lines show our model’s span of the /1-profile’s parameter space, compared to the results of MB15. Our full model calibrated to 
the simulations has pulled down the model’s exponential falloff, /3 from our earlier toy model, inferring an average gaseous halo 
model denser than the results of MB15. 


effect as a multiphase ISM, further motivating our listed 
truncation radius errors. 

Stellar Surface Density: - The GG model relies on the 
stellar surface density distribution, which we treat as 
an exponential disk (see Equation [l]) with characteris¬ 
tic ma ss M*, radial scale radius a+. TRoediger fc Henslerl 
(12005T) found the vertical structure and thickness of the 
disk did not impact the stripping dynamics, and thus we 
ignore uncertainty in &*, the scale height of the phys¬ 
ical density profile. The surface density is proportional 


to the stellar number density distribution. Ivan der Marell 
(120011) (Section 3.1) fit an exponential falloff to this num¬ 
ber density to find a scale radius ra = 1.44°. Starting 
with the data in their Figure 4 we used emcee to find a 
standard deviation in this value of ss .04°. This combined 
with the known 5% uncertainty in distance to the LMC 
leads to a combined error in our stellar scale radius of 6%, 
or a* = 1.7±.l kpc. Uncertainty in the total stellar mass 
is dominated by uncertainty in the stellar mass-to-light 
ratio, which vdM02 gives as M/I = 0.9 ± 0.2, and thus 
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a 22% error, or M* = 2.7 ± .6 x 10 9 Mq. We include an 
estimation of both these uncertainties by sampling from 
independent Gaussian profiles for both these variables. 
Sampling from these distributions holding all other vari¬ 
ables fixed allows us to gauge how sensitive our results 
are to these uncertainties. Using the observed trunca¬ 
tion radius of 6.2 kpc and the analytic Gunn / Gott 
prescription, we find the uncertainty in the scale radius 
introduces a 4% uncertainty in the implied pericentric 
halo density, n p , while the uncertainty in mass adds 23%. 
The latter is unsurprising, since in the Gunn/Gott model 
n p oc E* oc M*. 

Disk Orientation: — The LMC’s disk orientation (posi¬ 
tion and inclination angle) can potentially introduce er¬ 
rors on both the simulated and observed end of our con¬ 
siderations. Regarding the simulations, Section 14.51 ex¬ 
plored the role of disk-orientation in altering the LMC’s 
present-day HI truncation radius. Specifically, Figure 
1101 showed runs with intermediate angles between the 
disk-plane’s normal and the wind vector (neither espe¬ 
cially close to face-on nor edge-on wind) featured trun¬ 
cation radii R t that evolved in a nearly identical fash¬ 
ion for a given halo model. This suggests uncertainty 
in the LMC disk plane’s orientation and the direction of 
the LMC’s pericentric velocity vector inferred from this 
via backwards-integrated orbits ought to have a negligi¬ 
ble impact on the uncertainty in our model’s mapping 
between pericentric halo density and observed trunca¬ 
tion radius. Regarding the observed HI distribution, the 
disk plane’s orientation also enters into our calculations 
when we fit a constrained ellipse to the leading edge col¬ 
umn contour, to infer a truncation radius (specifically, 
it determines the ellipse’s axis ratio and orientation). 
vdM02 places the uncertainty in the position angle at 
10° whereas the inclination angle’s uncertainty is roughly 
half that. The spread in inferred truncation radius in¬ 
troduced by this is significantly smaller than uncertainty 
introduced by the LMC COM distance measurement, 
which is already in turn much smaller than the <r = .25 
kpc we introduced into R t ’s observed value due to the 
flocculcnt appearance of the leading edge. 

LMC Dark Matter Potential: — Our parameter choices 
for the DM halo were motivated by matching the ob¬ 
served LMC disk’s rotation curve while also matching 
our Burkert profile closely to an NFW profile with con¬ 
centration parameter c = 10. Our simulated rotation 
curve peak of 80 km/s is a dmittedl y on the l ower end 
of the observations (VdM14. lOlsen ct al.ll201llh An al¬ 
ternative DM potential impacts the dynamics in a few 
notable ways. First, a more (less) massive LMC would 
make dynamical friction more (less) efficient, altering its 
orbital behavior. Our sampling of LMC/MW masses 
described above accounts for this uncertainty. Second, 
an enhanced (diminished) circular velocity may make 
stripping more (less) efficient where the disk gas motion 
parallels the headwind (momentum-wise it has a “head 
start”). And third, the gravitational restoring force of 
the DM potential may alter the required ram pressure 
for stripping at a given radius, an issue increasingly rel¬ 
evant in dwarf systems. The GG model was first de¬ 
veloped to handle ~ 10 12 M 0 late-type systems falling 
into clusters. For these systems, it’s a good approx- 



Figure 15. A comparison of simulations with disparate LMC 
disk circular velocities (i.e. DM potentials), plotting quanti¬ 
ties as in Figure [9] Upward triangles, squares and downward 
triangles denote circular velocities of 100, 90 and 80 km/s, 
the final value representing our standard model throughout 
this work. Colors denote different MW gas halo core den¬ 
sity models. Across n p and wind models, we fail to find a 
significant trend with Vdic, and thus the DM potential does 
not appear to play a strong role in the setting the truncation 
radius. 

imation to assume a (thin) stellar disk dominates the 
restoring force during RPS (see Section [4761 for a discus¬ 
sion of the gas’ self gravity). For the LMC, the vertical 
acceleration of the DM potential at the truncation ra¬ 
dius r ~ 6 kpc is an appreciable fraction of the stellar 
surface density’s acceleration field, which suggests DM 
could alter the predicted truncation radius. These latter 
two issues have thus far not been captured in our models 
nor our simulations, and could introduce additional un¬ 
certainty. To explore whether or not this was the case, 
we ran ten additional simulations: five with a rotation 
curve peaked at 90 km/s and five at 100 km/s. We chose 
three n p € [0.8,1.2,1.7] x 10~ 4 cm~ 3 and at the median 
n p employed our three distinct halo core models. Fig¬ 
ure [T5] shows these results on our standard i?trunc vs. n p 
plot. For each of the five wind scenarios, the three dis¬ 
parate DM potential runs appear in a roughly random 
order, suggesting no strong trend with circular velocity 
(a trend would involve the symbols aligned as in the leg¬ 
end in the lower left-hand of the figure). Further, the 
scatter among these variations is less than the spread 
in observed truncation radius and that inherent to the 
simulations. Thus we ignore this uncertainty. 

Table 5.5 summarizes results of the sampled models, 
including a breakdown of spread in all the parameters 
considered. From these results, we find the uncertainty 
in the observed truncation radius of the leading edge ac¬ 
counts for the most scatter in our results: ~ 40% of 
the variance. We emphasize that our choice of variance 
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Inferred MW Halo Density, 

n p (48 kpc) 


Parameters Varied 

Mean 

Median 

Std Dev 

Simulation Scatter 

1.066 xl0 -4 cm -3 

1.044 

.2214 

Stellar Surf Density 

1.104 

1.081 

.2294 

Truncation Radius 

1.087 

1.047 

.2924 

Pericentric Velocity 

1.059 

1.053 

.1270 

Combined Errors 

1.118 

1.030 

.4950 


Table 6: Inferred MW gaseous halo density, n p at r p = 
48.2 ± 2.5 kpc from Sgr A*, from our GG model with an 
offset calibrated from our suite of wind-tunnel simulations. 
The top row incorporates only the scatter in the relation be¬ 
tween the HI truncation radius, Rt, and n p found from the 
simulations; the next includes only the mean offset from the 
simulations and spread in parameters for our modeled stel¬ 
lar disk; the next includes only spread in the observed Rt 
(assuming a Gaussian spread, although this choice does not 
strongly influence the results); the next row from including 
all various LMC pericentric orbital speeds from our 10,000 
orbit samples; and the last combining all these uncertainties. 

for this observed truncation radius was set to a conser¬ 
vatively large value of .25 kpc. Next, the uncertainty 
introduced by scatter in the simulated truncation radius 
(mostly due to oscillations of the disk edge) and impreci¬ 
sion in the stellar surface density profile each contribute 
roughly equally to the variance, ~ 25% each. Finally, 
the spread in pericentric velocity from our sampled or¬ 
bits contributes only ~ 8% to the uncertainty in our 
inferred value of n p . 

6. DISCUSSION 

6.1. MW Gaseous Halo Mass Estimate 

Baryonic mass accounts for a fraction = .17 

of matter in the universe. But stars and the ISM in 
the disks of L ~ L* galaxies account f or a far smaller 
share of their host halo’s total mass (e.g. iBehroozi et al.l 
HmJ. Although galaxies may efficiently disperse the re¬ 
maining baryons into the intergalactic medium (IGM), 
observations of that material have thus far fai led to ac¬ 
count for a sufficient mass in baryons dProchaska et al.l 
120111 ICen fe Ostrikedli999D . The MW’s CGM has of¬ 
ten been invoked as a reservoir of these so-called “miss¬ 
ing baryons”. MB13’s results suggests the MW hot 
halo accounts for 4.3 ± .9 x 10 lo Mm or roughly 50% 
of these baryons, though iGupta et alJ d2012i ) suggest a 
much larger mass of hot gas. 

Surveys of ~ L* galax ies at low-redshift sug gest a sim¬ 
ilar role for X-ray gas (lAnderson et al.l 120131) . The re¬ 
cent COS-halos results of QSO absorption in UV sug¬ 
gests the CGM of ~ L* systems is indeed rife with 
baryons, though at far lower temperatures: potentially 
half the galaxy’s mea n baryon budget may exist as diffuse 
T < 10 5 K CGM gas dWerk et al.ll20l4lTumlinson et all 
120131 ): while warm-hot 10 5 — 10 7 K mater ial could like¬ 
wise harbor a signifiant frac tion of mass (|Peenles et all 
120141 iTumlinson et al.|l201 lh . 

Our work thus far has focused on inferring a local den¬ 
sity measurement at the LMC’s last pericentric passage, 
and the analysis of Section [4751 showed this measurement 
is largely insensitive to the broader halo profile. In this 


section, we consider the range of total CGM mass esti¬ 
mates implied by this density measurement. Such an ex¬ 
ercise is indeed an extrapolation; in addition to invoking 
a specific density profile, we’ll be integrating well beyond 
LMC distances, assuming a consistent exponential falloff 
and spherical symmetry. 

For the density distribution we again assume a /3- 
profile. For each n p sampled as in Section [5J we also ran¬ 
domly sample a halo core density, nor , drawing from a 
uniform distribution of nor^P £ [.005, .003] correspond¬ 
ing to the range in MB 15 though with the upper bound 
set a factor of three higher, since their work fails to probe 
cooler phases of the CGM. We then infer /? from this core 
density and the sampled pericentric density. 

The left-hand panel of [Hi] shows the average density 
profile and 1 — a bound at each radius from these samples. 
The right-hand panel integrates this profile to find the 
total mass enclosed by each radius. From here, we find 
the standard runs imply a total CGM mass of 2.68 ± 
1.4 x 10 10 Mq. To place this bound in context, we can 
normalize to the “expected” mass of baryons in the MW, 
by multiplying the MW’s total halo mass by the mean 
cosmic baryon fraction, flb/f Im = -17. For a MW mass 
of 1.5 x 1O 12 M 0 (e.g. iBovlan-Kolchin et all 120131) . our 
standard runs imply the CGM contains ~ 10% of the 
Galaxy’s baryons. 

Table 6 quotes values at 300 kpc and also expresses 
these masses as a fraction of “available” baryons, i.e. 
.17Mhaio for a set of MW masses M £ [.7,1.0,1.5,2] x 
10 12 . From our work, the most reasonable estimate 
places the CGM baryon fraction at « 10% — 20% of the 
available baryons. In comparison, the COS-halos team 
finds anywhere from 25 — 45% in a cool, < 10 5 K form, 
and warm-hot contribution of 5 — 37%. MB15 in contrast 
predicts a < 50% contribution from X-ray gas. Our result 
is on the conservative end of these predictions, despite 
adopting wide bounds in the halo model and probing the 
full temperature span of gas. 

Although our work focused on a spherically symmet¬ 
ric /3-profile, our analysis has shown that the currently 
observed truncation radius is sensitive only to the den¬ 
sity at perigalacticon, where ram pressure is strongest, 
and thus our density bounds are consistent with a broad 
array of models for the diffuse CGM, so long as density 
tended to rise appreciably as the LM C moved deeper into 
the M W’s potential. Alternatively, iTonnesen fe Hr van! 
(I2008D examined RPS dynamics in a larger, forming clus¬ 
ter environment, where they found ram pressure experi¬ 
enced at a single radius from the cluster center can vary 
significantly. They found stripping to be a fast enough 
process in such an environment that truncation of the 
cluster galaxies was in fact a probe of dense pockets. If 
we assume our scaled down setting would obey similar 
dynamics, this implies our halo density measurement is 
in fact an upper limit, as the current LMC truncation 
radius could be set by the highest ram pressure pocket 
encountered by the cloud, rather than at perigalacticon. 

6.2. SMC Stripping Consistency Check 

The LMC’s companion galaxy, the Small Magellanic 
Cloud (SMC), provides another late-type dwarf travers¬ 
ing the CGM at high speeds at intermediate distances 
from the galactic plane. In this section we briefly verify 
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Figure 16. Density and enclosed mass profiles for the MW CGM, generated by sampling the diffuse CGM density at the 
LMC’s pericentric passage, from our stripping simulations, and then sampling /3-profile core densities from nor 8/3 £ [.005, .03]. 
Table 6 summarizes the mass enclosed and provides implied baryon fractions. 


Extrapolated Gas Halo Mass 


M gas (300 kpc) 

.7 x 10 12 

AW(.17 

1 x 10 12 

io) 

1.5 x 10 12 2 x 10 12 

2.68/^ x 1O 1U M 0 

21.58/^% 

15.11+^% 

10.07/^% 7.55l^% 


Table 6: The total CGM mass enclosed within 300 kpc from our LMC stripping calculations. Also shown are the CGM’s 
contribution to the Milky Way “Baryon Budget” defined as the halo’s total mass multiplied by the cosmic mean baryon 
fraction, = .17. From here we find the CGM likely contributes « 5% — 30% of the MW’s expected baryon mass. These 

results assume a spherical /3-profile calibrated to the CGM density found at LMC pericenter. 


Quantity 

Value 

Unit 

M* (2 kpc) 

3.1 x 10 8 

Mq 

a* 

.65 

kpc 

Mg as (2 kpc) 

4.0 x 10 8 

Mq 

&gas 

max, gas 

.65 

9.98 x 10 21 

cm -2 

ODM 

27.5 

km/s 

rsMC 

60 

kpc 

VSMC 

217 

km/s 


Table 7: SMC model parameters, as explained in Section[62] 

that our results are consistent with the survival of the 
SMC’s observed HI distribution. 

The SMC system’s position, orientation and motion 
are considerably less well-constrained than the LMC, 
rendering conclusions we arrive at here necessarily more 
speculative. In particular, the past pericentric approach 
of the SMC to the MW is highly uncertain as its orbit is 
strongly perturbed by the LMC. Furthermore, the sys¬ 
tem is surrounded by HI at appreciable column along 
all sight-lines, with no leading edge-analogue where RPS 
would be an obvious shepherd of HI (i.e. there is no loca¬ 
tion where gas is truncated relative to the stellar distri¬ 


bution as for the LMC). Thus RPS arguments are only 
useful in verifying that our density bounds allow for the 
survival of the observed SMC gaseous system. 

We employ the same analytic GG model as for the 
LMC. Table 7 summarizes all inputs for this model. We 
once again employ Equations 1 and 2 to describe ex¬ 
ponential profiles for the stellar and gaseous disks. We 
adopt a gaseous mass M gas (2 kpc) = 4 x 1O 8 M 0 and 
central surface brig htness of E maSigas = 9.98 x 10 21 cm~ 2 
(jBriins et al.ll2005l) . which for an exponential profile im¬ 
plies a scale radius a gas = .65 kpc. For the stars, we 
adopt a mass within the same radius of 3.1 x 10 s M 0 
(jStanimirovic et ahll2004l : Ivan der Marel et al.ll2009H . An 

accurate account of the gas’ restoring force must in¬ 
clude the SMC’s DM mass, which is relatively large 
and concentrated. We adopt a dispers ion velocity 
od m = 27.5 k m /s ([Harris fc Zaritskvll2006ll . and follow 
iGrcevich fc PutmanP 20091 1 bv modeling its affect via 

(14) 


UpUp ~ O [,Tl gas ( /f'fr 

unc) 


where the ambient gas’ ram pressure is formed from n p 
and v p , the CGM density and relative velocity at the 
SMC’s pericenteric passage. We examine the present 
properties of the SMC, which is moving with a Galacto- 
centric velocity of 217 ± 26 km/s (K13) at a distance of 
r smc ~ 60 kpc. For n gas we take the SMC ISM’s gas den¬ 
sity at the disk’s mid-plane (i.e. an upper bound). The 
final input is a truncation radius for the SMC’s gaseous 
extent, which we set to 3 .5 kpc, taken from the HI disk’s 
observed angular extent dBriins et al.ll2005l) . This analy- 
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sis implies a density upper bound just above 10 -4 cm -3 . 
This limit is consistent with the LMC fiducial case we’ve 
considered throughout this paper, but excludes the gas- 
rich scenario of Section 14.61 

6.3. Contributions to and Constraints from the 
Magellanic Stream 

The Magellanic Stream has itself been used to place 
constraints on the local halo gas density. Observations 
of the radius and velocity width of clumps of HI in the 
Stream can be used to derive their internal pressure. If 
the clumps are assumed to be entirely confined by gas 
pressure (no contribution from magnetic fields or turbu¬ 
lence), the necessary halo gas density can be estimated 
(jStanimirovic ct al.ll2002l : [Murall[2000). Assuming a halo 
temperature to 3.5 x 10 6 K. lAndcrson fc Brcgman (|201 Of ) 
place an upper limit on the hot halo electron density of 
n(45kpc) < 9 x 10 -5 cm -3 . Note, however, that the 
true distance to clumps i n the Stream is li k ely signifi¬ 
cantly larger than 45 kpc dBesla et al.ll2012f) . iHsu et ahl 
(1201111 conduct a similar analysis for clumps at the tip of 
the Magellanic Stream assuming a distance of 120 kpc. 
They find that pressure confinement arguments require 
halo densities from 10 -44 to 10 -4 ° cm 3 . Such results 
imply a substantial flattening of the density profile at 
radii beyond 50 kpc, inconsistent with a /3-profile. 

A pure ram pressure stripping origin for the Mag¬ 
ellanic Stre am has been ex p lored by many au¬ 
thors (e.g.. iMeu rer et al. 1 119851: iM oore fc Davis! 11994 
iHeller fc Rohlfsl 11994 iMastrouietro et al.l 120051 1 . Such 
models naturally explain the apparent lack of stars in the 
trailing Magellanic Stream, but do not explain the exis¬ 
tence of the Magellanic Bridge connecting the LMC and 
SMC or the existence of the Leading Arm (which leads 
the entire system). It is thus clear that tidal interac¬ 
tions between the LMC and SMC must play a dominant 
role in shaping the gaseous Magellanic System. However, 
this does not imply that ram pressure stripping provides 
a negligible contribution to the Stream. Indeed our anal¬ 
ysis illustrates that ram pressure will alter the structure 
of the LMC’s gaseous disk. In particular, while the bulk 
of the stream has a very l ow metallicity (0.1 Z 0 ), consis¬ 
tent with an SMC origin dFox et al.ll2013lh closer to the 
Clou ds the metallicity increases to 0.5 Z 0 (iRichtcr et al.l 
120131! . These observations suggests that the Stream is 
being enriched, possibly from gas recently ram pressure 
stripped from the LMC. 

We do not explicitly track the multiphase structure of 
the ISM in our simulations, instead we identify neutral 
HI by following only the evolution of gas in cells with 
densities above p = .03 cm -3 and column densities of 
10 19 cm -2 . Gas obeying these conditions is not found 
at large radii beyond the LMC disk, indicating that the 
LMC is not a dominant contributor to the Magellanic 
Stream. However, it is certainly possible that gas dis¬ 
obeying such conditions may in fact cool radiatively to 
form neutral HI in a full treatment of the multiphase 
ISM. As such, in Figure [T7] we plot the column of all 
gas that was once in the LMC disk in our fiducial simu¬ 
lation after 1 Gyr of evolution (regardless of density or 
temperature). The bottom panel illustrates the observed 
distribution of HI in the Magellanic System. We find 
that the bulk of the gas removed from the LMC exists at 
low column densities (< 10 19 cm -2 ). Furthermore, this 


material almost exclusively follows the past orbital tra¬ 
jectory oFthe LMC on th e plane of the sky. As pointed 
out bv lBesla et al.l (|2007f ) and K13, HST proper motions 
rule out scenarios where the past orbital trajectory of 
the LMC is coincident with the location of the Magel¬ 
lanic Stream on the plane of the sky. It is instead offset 
by ~ 10°. Such offsets are natural consequences of the 
tidal stripping of rotating disks; since MW tides are too 
weak to affect the LMC, this results favors a scenario in 
whic h the Stream forms thro ugh tidal stripping of the 
SMC IBesla et al.l (l20H l20Tnlf . 

The exact contribution of this low column density ma¬ 
terial to the mass budget of the Stream is speculative, 
as it depends sensitively on the original extent of the 
LMC’s gaseous disk and also the dark matter mass of 
the LMC, which controls how much gas can escape its 
potential well. For our fiducial simulation presented in 
Figure [T71 we find that ~ 7 x 1O 6 M 0 of gas (at any 
temperature or density) is removed from the LMC’s gas 
disk outside a radius of 9 kpc (the nominal extent of the 
LMC’s stellar disk). As such, we estimate that the LMC 
may contribute less than 1% to the total mass budget of 
gas outside the Magellanic Clouds (the total gas mass, 
including the ionized compone nt, is 2 x 10° Mm, a ssuming 
an average distance of 55 kpc IFox et al.l (i2014ll i. If all 
of this material were neutral HI, it would contribute no 
more than 6.4% of the mass budget of the head of the HI 
Magell anic Stream (1. 1 x 1 0 8 M(p, assuming a distance of 
55 kpc lPutman et al.ll2003H . 

Note that, in our fiducial model, the total LMC mass 
within a radius of 100 kpc is ~ 5 x 1O 1O M 0 . Given the 
observed gas and stellar mass of the LMC, a reasonable 
baryon fraction (3-5%) would imply a total halo mass 
of order 2 x 10 n M^ ( i.e., expectatio ns from abundance 
matching iBovlan-Kolchin et al.ll201 ll ). In such a massive 
halo it is likely that much of the stripped material would 
remain bound to the LMC. Our fiducial model is thus 
an optimistic estimate of the contribution of material 
ram pressure stripped from the LMC to the Magellanic 
Stream. Even in our lower LMC halo model, material 
that is removed from the LMC does not travel to large 
distances (the Stream spans more than 100° on the sky). 
As such, the vast majority of the Stream is not likely 
to be ram pressure stripped material from the LMC and 
this material would only be found close to the clouds 
consistent wi th the higher metallicity material fo und in 
observations (|Fox et al.ll2013l iRichter et al.ll2013| l. 

In all simulations we find that stripped gas is initially 
removed asymmetrically, building up in the North-West 
quadrant (upper right) of the LMC disk. This occurs 
because the clockwise sense of rotation of the LMC disk 
allows material to be removed more easily in the north, 
where disk gas is alre ady moving in the same dir ection 
as the headwind (e.g., iRoediger fc Bruggenll2006 1. This 
leads to a very specific prediction: gas in this region 
should have metallicities consistent with that of the outer 
disk of the LMC and higher than that of the tip of the 
Stream. This can be tested with QS O a bsorption line 
studie s such as that of lFox et al.1 ( 20158 ! and lLehner et al.l 
(|2012f l. Further observations to confirm the metallicity 
break in the Strea m observed with one QSO sightline 
near the Clouds by iRichter et al.l (I2013H are also war¬ 
ranted. 
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Figure 17. Upper Panel: gas column density distribution of the simulated LMC after 1 Gyr of evolution. No density or 
temperature cuts are applied. Lower Panel: the true distribution of HI gas in the Magellanic Clouds, Bridge and Stream in 
the same viewing perspective. Overplotted in grey scale is the distribution of stars in the LMC disk. The simulated material 
stripped from the LMC is of much lower column density (< 10 19 cm -2 ) than the bulk of the observed stream and follows the 
past orbit of the LMC on the plane of the sky, which does not trace the location of the Magellanic Stream. The exception is 
material in the North-West quadrant of the LMC (upper right of the LMC disk), which may be higher column. There is HI at 
high column observed at that location, which we predict should have originated exclusively from the outer disk of the LMC. 
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7. SUMMARY 

We have provided a density measurement of the Milky 
Way’s circumgalactic medium at r = 48 kpc by exploring 
the dynamics of ram pressure stripping of the Large Mag¬ 
ellanic Cloud via analytic arguments and a large suite 
of three-dimensional hydrodynamic simulations. RPS 
probes the diffuse CGM in all its phases, providing a 
robust measurement independent of line studies in emis¬ 
sion and absorption. We adopted a simple model for 
the LMC’ s stars, gas an d DM profiles and adopted the 
model of iGunn fc Gotti (11972H (GG model) to explore 
how a ram pressure headwind truncates the HI profile of 
disk systems, providing a wind computed from models 
for the LMC’s orbit through the MW potential and sam¬ 
pling from a /3-profile for the CGM’s density. Our main 
conclusions are: 

1. The presence of stars well beyond the HI truncation 
radius of the Large Magellanic Cloud’s leading edge 
suggests ram pressure stripping played a strong role 
during the LMC’s recent pericentric passage trun¬ 
cating its gaseous disk. 

2. From analytic arguments, assuming the LMC’s 
disk was truncated by the peak ram pressure at 
pericentric passage, we found the GG model pro¬ 
vides a direct map from the LMC’s observed trun¬ 
cation radius along its leading edge to a CGM den¬ 
sity at r ss 48 kpc from the Galactic Center. This 
toy model implied a halo density « 5x 10 -5 cm -3 . 

3. From our suite of simulations sampling a broad 
range of halo models, we found the LMC’s lead¬ 
ing edge truncation radius was indeed tightly cor¬ 
related to the MW CGM density at pericentric pas¬ 
sage, as predicted in the GG model, and roughly 
independent of the following considerations: 

(a) The broader halo profile’s core density and 
power-law falloff - i.e. the density measure¬ 
ment was very localized to pericentric passage, 
thanks to the steep rise in ram pressure at 
this point. This implies our measurement is 
independent of our choice to use a /3-profile 
and applies to any distribution for the diffuse 
CGM. 

(b) The disk orientation with respect to the wind, 
provided it was neither head-on nor precisely 
face-on, making our result insensitive to the 
precise orientation of the LMC disk at peri¬ 
centric passage. 

(c) The DM potential and rotation curve peak of 
the LMC 

(d) The restoring force of the disk in such a gas- 
rich late-type dwarf is best described by in¬ 
cluding not only the stellar disk’s gravita¬ 
tional potential but also the self-gravity of the 
gas. 

4. The simulations also elucidated some disagreement 
from/extensions of the GG model, including: 

(a) The presence of extra-planar material in the 
stripping zone renders the truncation radius 


seen in column viewer-dependent, as mate¬ 
rial in the process of liberation from the disk 
contributes gas at high column depending on 
where one views the galaxy from. This effect 
can introduce a substantial asymmetry in the 
observed extent of the disk in different quad¬ 
rants of the system. 

(b) Our most relevant simulations included a sig¬ 
nificant 5% offset from the GG model, i.e. our 
truncation radius for a given headwind ram 
pressure was 5% larger than the analytic pre¬ 
diction with the gas’ self-gravity. A portion 
of this offset was present regardless of viewing 
perspective (from the solar line-of-site, face-on 
and edge-on projections of the disk), and at¬ 
tributable to time-varying oscillation’s in the 
disk’s extent. 

5. Analysis of this suite of simulations, for our fiducial 
LMC disk model, yielded an implied CGM density 
substantially higher than the analytic toy model, 
with n(48.2 ± 2.5 kpc) = 1.1 ± .45 x 10 -4 cm -3 . 
The error bounds on this include a host of modeled 
uncertainties, including those related to the scatter 
in the simulations (mostly due to the oscillatory be¬ 
havior of the disk edge), orbital uncertainties (in¬ 
cluding that of the LMC/MW mass ratio), the floc- 
culent appearance of the LMC disk’s leading edge 
(the dominant source of error) and uncertainty in 
the stellar profile. 

6. From our simulated results, sampling parameters of 
a /3-profile and integrating to 300 kpc, our standard 
LMC model results yielded a total CGM gas mass 
of 2.68 ± 1.4 x 10 10 Mq, which could account for 
roughly 15% of a 1O 12 M 0 MW’s expected baryon 
content. 

7. It is unlikely that material ram pressure stripped 
from the LMC can account for more than a few 
percent of the mass budget of the Magellanic 
Stream. The simulations suggest that material to 
the North-West of the LMC’s disk consists of gas 
stripped from the LMC (but not necessarily un¬ 
bound). This material should have metallicities 
consistent with the outskirts of the LMC’s disk, a 
prediction that is testable through absorption line 
studies towards background QSOs. 

Movies for our fiducial simulation can be found at 
https: //lavinia.as.arizona.edu/ gbesla/MWHotHalo.html 
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APPENDIX 

LMC DISK PLANE VELOCITY CORRECTIONS 

To produce the HI velocity map in Section 14.5.51 we needed to subtract motion of the gas due to the LMC’s COM 
motio n and its solid body rotation. To do so, we employ the machinery laid out in vdM02 and Ivan der Marel fe Cioiiil 
(j2001t ). This appendix briefly describes the necessary equations along with a brief motivation of our parameter choices. 

This prescription takes as inputs the LMC’s COM position, (a o, <5o, .Do); disk orientation (*, 0); COM proper motion 
in the LOS frame, v x ,v y ,v z , (see Section l4~2ll : and the change in the disk plane’s inclination angle, di/dt. Note the 
change in position angle dO/dt does not contribute to the radial velocity. Values for all these quantities were taken 
from the aforementioned sources, and are summarized, for convenience, in Tabic [5] With these values in place, a set of 
angular coordinates, p and <j>, can be computed at every point in the projected image in terms of celestial coordinates: 

cos p = cos 5 cos S o cos (a — ao) + sin i5 sin 5q 
sin p cos <f> = — cos <5 sin(a — ao) 

sinpsin</) = sin(5 cos<5o — cos i5sin cos(a — ao) . (Al) 

although these are three equations in two unknowns, the system is not overdetermined since trigonometric functions 
are not one-to-one. From there, we can recast the LMC’s cartesian velocity vector into an angular form, still in the 
LOS frame: 



-v z \ 

0 * + ^) 1/2 
tan ~ 1 (v y /v x ) — 90°/ 


(A2) 


Here v sys is motion away form the observer, vt is speed of the motion in the plane of the image (transverse), and 0 t 
is the position angle of this transverse motion, measured counterclockwise from the y north-vector. Finally we’re in a 
position to compute the radial velocity correction, which must be subtracted from the observed velocities: 

uios(p, $) = v sys cos p + vt sinpcos(4> — ©t) + Do(di/dt) sinpsin($ — ©) , (A3) 

where 0 = 9 — 90 and d* = <f> — 90. The first term in this expression corresponds to the COM proper motion along 
the line of site, and differers from a naive subtraction of the bulk motion by the cosine factor. The second term is 
more pernicious, introducing a velocity gradient across the image plane of amplitude ~ 100 km/s, which applies a 
systematic twist to the observed velocity field. The final term is related to any rotation of the LMC disk plane about 
the line of nodes, which all evidence at present has failed to detect at a significant level (Ivan der Marel fe Kallivavalill 












